SCATTERING, HOMOGENIZATION AND INTERFACE EFFECTS 
FOR OSCILLATORY POTENTIALS WITH STRONG 
SINGULARITIES 



VINCENT DUCHENE* AND MICHAEL I. WEINSTEINt 

Abstract. Wc study one-dimensional scattering for a decaying potential with rapid periodic 
oscillations and strong localized singularities. In particular, we consider the Schrodinger equation 

H.4,= {-dl + V {x) + q{x,xfe))i> = k 2 il> 

for k £ R and e <C 1. Here, q(-, y + 1) = q(-, y), has mean zero and |Vb(^) + <?(^i -)| — > as \x\ — > oo. 
The distorted plane waves of H e are solutions of the form: ey±{x] k) = e^ lkx + Uj_(x;fe), «|_ 
outgoing as \x\ — > oo. We derive their e small asymptotic behavior, from which the asymptotic 
behavior of scattering quantities such as the transmission coefficient, t e (k), follow. 

Let t^ om (k) denote the homogenized transmission coefficient associated with the average poten- 
tial Vq. If the potential is smooth, then classical homogenization theory gives asymptotic expansions 
of, for example, distorted plane waves, and transmission and reflection coefficients. Singularities of 
Vb or discontinuities of g e are "interfaces" across which a solution must satisfy interface conditions 
(continuity or jump conditions). To satisfy these conditions it is necessary to introduce interface 
correctors, which are highly oscillatory in e. 

Our theory admits potentials which have discontinuities in the microstructure, q c (x) as well as 
strong singularities in the background potential, Vo(x). A consequence of our main results is that 
t c (k) — tQ° m (k), the error in the homogenized transmission coefficient is (i) 0(e 2 ) if q E is continuous 
and (ii) 0(e) if q t has discontinuities. Moreover, in the discontinuous case the correctors are highly 
oscillatory in e, i.e. ~ exp(2-7ri^-), for e -C 1. Thus a first order corrector is not well-defined since 
e- 1 (t e (k) - t$ om (k)) does not have a limit as e -> 0. This expression may have limits which depend 
on the particular sequence through which e tends to zero. The analysis is based on a (pre-conditioned) 
Lippman-Schwingcr equation, introduced in [9]. 
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1. Introduction. An important method for computing the effective properties 
of highly oscillatory media is the method of homogenization. The goal of homoge- 
nization is to approximate a highly oscillatory medium, described by a differential 
equation with oscillatory coefficients, by an approximate and homogeneous medium, 
described by a "homogenized" differential equation with constant or slowly varying 
coefficients. In its regime of validity, the homogenized differential equation (i) predicts 
effective properties which are approximately those of the heterogeneous medium and 
(ii) is, by comparison with the full problem, much simpler to study either analytically 
or by numerical simulation. 

While the homogenized limit can often be obtained by a formal multiple scale 
expansion or by variational methods [3l [101 Q] HI], these expansions are typically 
valid in the bulk medium, away from boundaries, discontinuities or more singular sets 
of coefficients. Indeed, solutions to elliptic operators with oscillatory coefficients on 
bounded domains have been shown to require boundary layer correctors, which are 
sensitive to the manner in which the microstructure meets a boundary [131 111! [2l [6j [7] 
or interface |15j . Furthermore, the importance of correctors to homogenization due 
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to interface effects, boundary layers etc. is explored analytically and computationally, 
in the context of accurate estimation of scattering resonances in [51 H] . 

In this article we study the scattering problem for the one-dimensional time- 
independent Schrodingcr equation 



The potential, V f ~(x) = Vo(x) + q(x,x/e), is the sum of a slowly varying part with 
smooth and singular components, Vq = V reg + V s i ng , and a rapidly oscillatory part, 
q e (%) = q (x, x/e) , e <C 1. V e {x) is assumed to decay to zero as x tends to infinity. We 
also assume V e (x) > 0, a simple way to restrict to the case where H e has no discrete 
eigenvalues (bound states) and has only continuous spectrum (extended / radiation 
states). The wave number, k, is fixed and we study the e— small behavior. 

Many physically important scattering properties are not captured by leading or- 
der homogenization. Line-widths and imaginary parts of scattering resonances are key 
to quantifying the lifetimes of metastable states in quantum systems, or in electro- 
magnetics, the leakage rates of energy from photonic structures; see [HI [H] and ref- 
erences therein. In jSJ [S] it was shown that inclusion of even the first non-trivial 
correction due to microstructure can yield large improvements in the approximation 
of such scattering quantities. Since, as we shall see, defects and singularities can be 
responsible for the dominant correctors and these contributions arc not captured in 
smooth homogenization setting, we therefore seek a better understanding of homog- 
enization for wave / scattering problems in their presence. In this paper we ask: 

How are scattering properties, such as transmission and reflection coefficients, t t (k) 
and r e (k), influenced by interfaces, defects and singularities? 

The heart of the matter is an asymptotic study of the distorted plane waves, 
solutions of (H e — k 2 )ip = of the form: 

ey±{x; k) = e ±lkx + u±(x; k), u s ± outgoing as \x\ — > oo, for e small. 

Consequences of our analysis include the following: 

1. Theorem 15.11 provides a convergent expansion of the distorted plane waves of 
Hq = — d^ + Vo + Q, which is valid for a large class of perturbing potentials, Q, 
which may be pointwisc large, but highly oscillatory (supported at high fre- 
quencies although not necessarily periodic). Thcorcm l5.5l is the corresponding 
expansion for the transmission coefficient t[k; Q]. By Proposition 15 . 31 we can 
apply Theorems 15.11 and 15 .51 to Q(x) — q e (x) = q (x,x/e), where q(x,y) is 1— 
periodic in y, decaying as |x| — > oo, and satisfies Hypotheses (V). 

2. Theorem 12.11 implies that: 



(i) t e (k) - 4 om (k) = 0(e 2 ) if q e is continuous and 

(ii) t e (k) - t% orn (k) = 0(e) if q e has discontinuities. 

For q e discontinuous interface correctors, which are highly oscillatory in e, 
enter the expansion; see the discussion in section [4] concerning failure and 
restoration of interface conditions at singularities of Vq or discontinuities of 
q e . These correctors are related to the asymptotics of boundary layers arising 
in work on homogenization of divergence form operators on bounded do- 
mains [131 HU ID H H]- Since these correctors involve e dependence of the 
form: - exp(27ri|), e < 1, ^ v 6 M, the expression e" 1 (t e (k) ~ t% om (k)) 




(1.1) 
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does not have a limit as e — >• 0, and a correction to the value of £g orn (/c) is not 
well-defined. However, there can be limits which depend on the particular 
sequences through which e tends to zero. See the more detailed discussion 
after the statement of Theorem 12.21 
Outline of paper: In section [2] we state detailed hypotheses and our main theorems 
on transmission coefficients, Theorems 12.11 and 12.21 which depend on our analysis 
of distorted plane waves (Theorem 15. ip . We also present the results of numerical 
simulations designed to illustrate the relationship between regularity of the potential, 
V e , and e small asymptotics of the transmission coefficient, stated in Theorem l2.ll In 
section [3] we present the technical background on one-dimensional scattering theory. 
In section [4] we derive, by including interface correctors to an expansion derived by 
the classical method of multiple scales, an expansion of the distorted plane waves and 
of the transmission coefficient valid to all orders in the small parameter e . Section [5] 
contains rigorous proofs of the expansion of the distorted plane waves (Theorem 15. ip 
and transmission coefficients (Theorem 15. II and 15. 5p with error bounds. The proof is 
based on the reformulation of the scattering problem as a pre-conditioned Lippman- 
Schwinger equation, an approach introduced in [3] . Appendix [A] contains a brief 
discussion of the numerical methods used in the simulations. Appendix [C] contains 
the technical proof of operator bounds which are central to the proofs in section [5] 
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fruitful discussions. VD was supported, in part, by Agencc Nationale de la Recherche 
Grant ANR-08-BLAN-0301-01. MIW was supported in part by NSF grant DMS-07- 
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the Courant Institute of Mathematical Sciences, where he was on sabbatical during 
the preparation of this article. VD would like to thank the Department of Applied 
Physics and Applied Mathematics (APAM) at Columbia University for its hospitality 
during the Spring of 2008 when this work was initiated. 

2. Main results and Discussion. We begin with the key hypotheses. Hy- 
potheses (V) make precise the decomposition of the potential, V, into regular, sin- 
gular and oscillatory parts. Hypothesis (G) specifies, for the cases of generic and 
non-generic potentials, Vq, the admissible values of the wave number, k. We then 
state and discuss our main results concerning the transmission coefficients, in the 
small e limit. 

Hypotheses (V) 



V e (x) 



Vq{x) + q e (x), real — valued 
(x) + q e {x), 

q(x,-), V'(x)>0, 



(2.1) 




(2.2) 



where 



1 . Singular part of V e , V s 



sing ■ 



N-l 




(2.3) 
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2. Regular part ofV e : V reg G L 1 ' 2 (R) with 



\\V\\ L i.2 = / (l + |s|) 2 \V(s)\ ds < oo. 



(2.4) 



3. Rapidly varying part of V e , q e (x) = q (x, 2 ): The mapping (x,y) H> q(x,y) 
is 

(a) 1— periodic, i.e. for each x G R, q(x,y + 1) = (7(2/), 

(b) mean zero with respect to y, i.e. for each x G R, 



9(35,2/) dy = 0, 



(2.5) 



(c) q G pC^.L 2 per , the set of functions g : R x S* 1 — > R, such that there 
exists a finite partition of R 



OO = Qq < di < (22 < • • • < dM < om+i = +00 



with 



XI / ll5(-,J/)llc3 ( „ dy < 00 
4. We shall work with the Fourier expansion of q(x,y), written as 
q(x,y) = £ qj (x)e 2 ^y, q 3 {x) = f e" 2 ™" <q(x,y) dy 

and assume 

I [ \q(x,y)\ 2 dy dx= V f \ qj {x)\ 2 < 00, 
Jr Jo 1 1^, Jr 

f \q(x,y)\ 2 dy = kj(x)\ 2 ">0, 



as — > 00 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



5. Proposition 15. 31 which is a step in proving Theorem 12. 11 requires more decay 
at infinity for q e : there exists p > 8 such that 

(l + l-ry/^GL 2 , |j|>l, and X;|(l + |.| a )^| <oo, (2.10) 

lil>i 



■j- ((1 + \x\ 2 yl 2 q {x)) G L 2 and sup 



dx 



b1>i 



< 00. 



L- 



(2.11) 



Hypothesis (G) If Vq is generic (see Definition 13. 6[) , then the wave number, 
fc G K, an arbitrary compact subset of R. If Vo is not generic, then the compact set 
K must be such that ^ K . 

Remark 2.1. J/ Vb is no/ generic (as for example Vq = 0), then the expansions 
we present in Theorem \2.1\ and Theorem \2.2[ are not uniform in a neighborhood of 
k = 0. This will be the subject of a future paper. 
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The aim of this article is to understand the scattering properties for this class of 
potentials. In particular, we are interested in the influence of combined microstruc- 
ture (q e ) and singularities (V s ,„j) on the reflection and transmission coefficients and 
distorted plane waves (see below). Formal application of classical homogenization 
theory (see for example [3]) suggests that the leading order (in e — >• 0) scattering 
behavior is governed by the averaged (homogenized) operator —d% + Vo(x); see (|2.ip . 
For example, if V e (x) is smooth (in particular, V s % n g = 0), then the transmission 
coefficient satisfies the expansion 

t e (k) ~ t% om (k) + et\ om {k) + eH 2 lom (k) + ... (2.12) 

where t^ om are computed from the formal 2-scale homogenization expansion. In par- 
ticular, t^ om is the transmission coefficient associated with the averaged potential 
Vq(x), However, homogenization is a theory valid only in the bulk, away from bound- 
aries or non-smooth points of coefficients. For our class of potentials, this expansion 
must be corrected. 

Our main result is the small e characterization of the distorted plane waves pre- 
sented in Theorem 15. II A key consequence of our analysis is the following: 

Theorem 2.1. LetV c (x) = V (x) + q e (x) withV and q € (x) = q(x,x/e) sat- 
isfying Hypotheses (V), and k G K a compact subset o/ffi satisfying Hypothesis (G). 
Denote by ey ±(x; k) the distorted plane waves associated with the unperturbed oper- 
ator —d 2 + Vq(x); see section^ 

Then, there exists eo = £o{K), such that for < e < eo, the transmission coef- 
ficient t e = t e {k) (see (|3.7|) ) associated with V e (x) satisfies the following expansion 
uniformly in k G K : 

f(k) = C m (k) + et\(k) + e 2 (#"»(*) + tl(k) ) + f rem (k), (2.13) 

where t^ om {k) denotes the transmission coefficient, associated with the average (ho- 
mogenized) potential Vq and 

*iO) = ^X]ev +(aj;fe)ev -(aj;fc) Ma, J— » ( 2 - 14 ) 

J =1 |i|>i 

t h 2° m (k) = ^Y1 1~ 2 I l*( z )! 2 e v a -(z; k)e Vo+ (z; k)dz, (2.15) 



a 2inl- 



tl(k)= TTT^Yl t d x( e Vo+fak) e Vo -{x;k)qi(x)) ] a . — ^ — , (2.16) 

tr em (k) = o(e 2+ ) 7 more precisely quantified in Proposition \5.6\ (2-17) 

(a) tj° m , j = 0,2..., denote the expansion coefficients for the transmission 
coefficient obtained from the two-scale (bulk) homogenization expansion, valid 
for smooth potentials. 

(b) t\ arises due to discontinuities in x i— > q(x, ■), and 

(c) t\ arises due to both the singular part of the potential, V s i ng , and discontinu- 
ities in x i-> q(x, •) or x H> d x q(x, •). 

t\ and t\ are uniformly bounded, for e small. However each is a sum over rapidly 
oscillating (as e — >• 0,) terms of the form exp(i^), corresponding to discontinuity points 
of q e , respectively points in the support ofV s i ng . 
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Theorem 12.11 is a consequence of the more general Theorem 12.21 stated below, 
which follows from the asymptotic study of the convergent expansion of the dis- 
torted plane waves, presented in Theorem 15.11 The proof of Theorem 15.11 is based 
on construction and asymptotic study of the scattering problem via a pre-conditioned 
Lippman-Schwingcr equation. This approach is quite general and applies to the per- 
turbation theory of Schrodinger operators of the form 

H = -d 2 x + V (x) + Q(x), 

where Q is small in the sense that HIQIII - (7- A)-tQ(7- A)"* is small. 

This formulation was introduced in [9] to study the perturbation of scattering res- 
onances due to high contrast microstructure perturbations of a potential. If Q is 
a "microstructure" , roughly meaning that it is supported at high frequencies, then 
HIQIII is small. Here, we apply this method and obtain a convergent expansion of 
Q i— > ev +Q(x, k) for fixed k and |||Q||| sufficiently small. The expansion of the trans- 
mission coefficient, Q i— > t\/ a +Q(k), is a direct consequence of: 

Theorem 2.2. Let V{x) = V (x) + Q(x) with V satisfying Hypotheses (V) 
and (1 + \x\ 2 ) p / 2 Q £ L 2 , for p > 8. We use the following norm on Q, see section [5.2[ 



(Do)- 1 (1 + M 2 ) P/4 Q (l + \x\ 2 ) p/4 (D Q y 

Set k £ K a compact subset o/R satisfying Hypothesis (G), and denote by ev ±{x; k) 
the distorted plane waves associated with the unperturbed operator — d 2 + Vo(x); see 
section^ Denote by t — t(k,Q) = t(k) the transmission coefficient (see (|3 . T[) ) asso- 
ciated with V(x). There exists To = to(K) such that for < |||<5||| < To(K), we have 
the following expansion which holds uniformly in k £ K : 

t(k,Q) = t' lom (k) + hlQ] + t 2 [Q,Q] + t rem (k), (2.18) 

with tQ° m (k) the transmission coefficient, associated with the average (homogenized) 
potential Vo, and the following: 

tl[Q] = 27ifc/ Q(0 ev +(C;k) e Vo -((;k) d(, (2.19) 

h[Q ' Q] ^^kJ Q Rv (k)(Q(() e Vo+ ((;k)) e Vo _(C,k) d(, (2.20) 
trem(k) = O (| | |Q| 1 1 2+ ) and more precisely estimated in Theorem \5.5\ (2.21) 

Here, Rv {k), ev +(x;k) and ey _(x;fc) being defined in section^ 

Remark 2.2. Symmetry considerations: There is a class of potentials, q e , 
whose members are discontinuous, and yet the ( oscillatory in e) correctors, tj , j > 1 
vanish. In subsection \2.1\ we explore families of such structures. Indeed, let us apply 
Theorem \2.2\ with V = V e satisfies Hypotheses (V) as well as the additional properties: 
Vo even and q e "separable" : 

Vo(x) = V (-x), q e (x) = q {x) q per . 

One can easily see that Vo is even implies that ey +(-; k)ev -(-', k) is even. Therefore, 
if qo and q per are of opposite parity, then x i— > ev +{x; k) ey n _(cc; k) q e (x) is odd and 



SCATTERING FOR OSCILLATORY POTENTIALS WITH STRONG SINGULARITIES 7 

therefore ti [q e ] (fc) = for any e > 0. It follows that for such potentials, and even if q e 
is discontinuous, the leading order correction to tQ° m (k) is t 2 [Q,Q] which is of order 
C(e 2 ) (see section \5.4\ l- Moreover, in this special case the second order corrector is 
well defined: 

lim e- 2 ( t e (k) - 4 om {k) ) = 4 om (k). 



The three subplots of figure 12.11 illustrate the results of Theorem 12.11 on the be- 
havior of t e — to for several contrasting choices of potential V e = Vb + q e , where Vb is 
a finite sum of Dirac delta functions, at equally spaced points^ 

• The left panel of figure [2~T1 corresponds to the case where q e is discontinuous. 
It shows that 

t e _ t hom = Q ^ e _^ and yet e -l ( t e _ f ho m ) doeg nQt haye & jj^ 

• The center panel of figure 12.11 corresponds to the case where q e is a smooth 
function, and Vq is a Dirac delta function. Here, 

f - 4 om = C(e 2 ), e -> 0, and yet e" 2 ( ^ - t^ om ) docs not have a limit. 

• The right panel of figure 12.11 corresponds to the case where q e is a smooth 
function, and Vb is a smoothed out Dirac delta function. Here we find 

t e _ t hom = 0(e 2 ), e^O, and \ime~ 2 ( t £ - t' lom ) is well - defined. 

This phenomenon of indeterminacy of higher order correctors, due to boundary layer 
effects is discussed, in the context of a Dirichlet spectral problem [T51 111) . 




Fig. 2.1. Illustration of Theorem \ 2. H via plot of log \t e — t^° m \ versus loge -1 for the case of 
q discontinuous and Vb a sum of Dirac delta Junctions (left panel, average slope 1), q smooth and 
Vb i sum of Dirac delta functions (center panel, average slope 2 ). The right panel (slope 2) is for 
the case where V e = Vb is a smooth approximation of a finite sum of Dirac delta- functions. 

The transition between the cases of a regular potential and a potential containing 
singularities is illustrated in Figure |2~2"1 The three panels show the behavior of t e — 1 
with respect to e, where the potential V e = Vo + q e satisfies q e is smooth and Vq is 
a sum of smoothed out Dirac delta functions. From right to left, Vq is an improving 
approximation of Dirac delta functions. 



The precise functions and parameters used to obtain the plots displayed in figurcs l2,ll and [2.2l 
are given in Appendix [^] page [35] 




2.1. Some specific structures. We now study in detail two natural and illus- 
trative classes of potentials: 

1. We first consider a one-parameter family of structures, which are truncations 
of a smooth potential, where for certain parameter ranges the manner of trun- 
cation causes a discontinuity. The latter corresponds to cleaving a periodic 
structure in a manner not commensurate with the background medium: 

= cosf^+0^ l[_i,i](a?). (2.22) 

We arc obviously in the case related in Remark 12. 2[ with Vq = (so that 
ey±(x;k) = e ±lkx and t^ " 1 = 1). More precisely, it is easy to show that 

tl(k;6) = — ^2e Vo+ (a j ;k)e Vo -(a j ;k) ^ [qi] aj — - — 
n i=i |J|>i 

= _ cos(e) sin^-j. 

In general, t\(k) ^ but for 9 = ^+mn, raeZ, q e (-; 9) is even and therefore 
for all k £ M. and e > 0, we have ti[q e ](k) — 0. 

2. Our second example is a piecewise constant (discontinuous) structure which 
is smoothly truncated 

V^x-9) ee h per {^+9) e-c-iH«+D (x), (2.23) 

with h per (y) the 1-periodic function such that h(y) = —1 for y £ (—1/2, 1/2], 
and h(y) = 1 for y £ (1/2,3/2]. 

Since the slow- varying part of q e (x) is smooth, and Vq has no singularity, 
Theorem 12.11 predicts that 

t e_ t hom = ^ e ^ Q) and lime- 2 ( f - t* om ) = t\ om is well - defined, 

even though the function q e (x) has internal discontinuities. In Figure [2~U1 we 
plot log \ t e — ig om | versus loge -1 for the two potentials V{ and V^, setting 
k = 1, and 9 = 0. 
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Fig. 2.3. Plot o/log|i e -t% om \ versusloge- 1 for the case of the potentials V{(x;8) in ll2~22l) 
(left panel, slope 1), and V^ix; 8) in 1 12. 2311 (right panel, slope 2). One has k = 1, and 8 = 0. 

3. Background on one-dimensional scattering theory. For simplicity, we 
consider potentials, W, which have no localized eigenstates, i.e. the spectrum of 
—d 2 + W(x) is continuous. We further assume that W has the form 

W = W reg + W sing , with 

W reg G L^+iR), 

N-l 

W S i n g = /J cj 5{x — Xj), where Cj,Xj € R, xj < Xj+i. 
j=o 

We now introduce an appropriate notion of solution to the Schrodingcr equation 

(H w ~ k 2 )u = (^--^ + W(x)-k 2 ^u = 0. (3.1) 

Let [£/]f denote the jump in U at the point £, i.e. 

[U} 6 = Um U{x) - lim U(x). (3.2) 

x— ►£+ a:— 

Definition 3.1. FFe say i/iaf u is a solution of time-independent Schrodinger 
equation (|3.1[) i/ u zs piecewise C 2 , and satisfies the (|3.ip on 1\ suppWgjng = K \ 
{xo, . . . , xat_i} as weH as i/ie jump conditions 

■ [u] x = 0, xeR, 

< =0 i/x eR\suppW si „ 5 , (3.3) 

[^w]^ = Cju{xj) where x 3 £ supp PV" S ;„ 9 . 

Of special interest are the Jost solutions, defined below. 

Definition 3.2. The Jost solutions f±(x;k) = m± (x; k)e ±lkx are the unique 
solutions of (|3.1j) . such that 

lim m±(x;k) = 1. 

x— f ±oo 
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This definition is valid, as we see in Appendix [Bj We shall use some smoothness and 
decay properties of these solutions, that are also postponed to Appendix [Bj for the 
sake of readability. 

With the help of the Jost solutions, we are able to define scattering quantities, as 
the transmission and reflection coefficients, and the distorted plane waves. 

Since f±(x; k) and f±(x; —k) are solutions of ()3.1[) . and are independent for k =/= 0, 
there exists unique functions t±{k) and r±(k), such that 

/-(*,*) = + ^ 

/+(*,*) = + /-(*,-*)■ 

It is then easy to check that t + (k) = t_{k) = t(k), and that t{k) and r±(k) are 
continous at k = 0. The distorted plane waves ew± (x; k) are then defined by: 

Definition 3.3. Given a potential W(x), we define ew±(x;k), the distorted 
plane waves associated with Hy/ by 

e + (x;k) = t(k)f + (x;k) = t(k)m + (x; k)e ikx (3.4) 
e-(x;k) = t(k)f-{x;k) = t(fe)m_(a; k)e~ ikx . (3.5) 



The distorted plane waves ew±{x; k) play the role for H\y that the plane waves 
e ±ikx p^ a y £ or _ —Q2^ as we see below. Let us first introduce the notion of outgoing 
radiation as \x\ — > oo. 

Definition 3.4. U(x) is said to satisfy an outgoing radiation condition or to be 
outgoing as \ x \ —> oo if 

( d x =F ik )U 0, as x — > ±co. 



Proposition 3.5. Given a potential W (x) , e\\r±{x]k), the distorted plane waves 
ew±{x;k) are the unique solutions of (|3.ip satisfying 

e w± (x;k) = e ±lkx + outgoing(x). (3.6) 

More precisely, they satisfy the following asymptotic relations [3]: 



e w+ (x;k) — (e lkx + r + (k)e~' lkx ) — > as x — >• — oo, 

ew+(x;k) — t{k)e lkx — > as x — s- +oo, 

ew-(x;k) — t(k)e~ %kx — > as x —> — oo, 

e W -(x;k) - (e~ lkx + r_(k)e lkx ) — > as x -> +oo. 



(3.7) 



A consequence of the relations (|3 . 7[) is the Wronskian identity: 
Wron (ew+(-; k), ew -{■'■, k)) = e w+ d x e W - ~ d x e w +ew- = -2ik t(k). (3.8) 
In terms of the Jost solutions: 

Wron (/+(•; k), /_(•; k)) = k^O. (3.9) 

t(k) 
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By analyticity in W, potentials for which Wron (/+(•; k), /_ (•; k))\ k=0 = are 
isolated in the space of potentials. 

Definition 3.6. A potential W is said to be generic if 

Wron(f + (x]0)J_(x]0)) = Wron (m+(x; 0), m_ (x; 0)) ^ 0. 

Otherwise, the operator Hyy is said to have a zero-energy resonance, i.e. H^u = 
has a non-trivial solution that is bounded both as x — > oo and as x — oo. 

Note that the potential W(x) = is not generic since m+(x; k) = rri-(x; k) = 1. 
If W is generic, we have [4| [T2 | fT5 ] 



2? k 

t(k) = . nw , m + o(k) = O(k), \k\^0. (3.10) 

Wron (/+(•; 0), f-(x; 0)) 

In particular, f(0) = and r±(0) = — 1. 

A simple calculation then yields the following expressions for the outgoing Green's 
function (resolvent kernel) and the outgoing resolvent, Rw(k), k ^ 0: 

{ -2ikt(k) e w-(y;k) e w+ (x;k), y < x, 
(3.11) 
-2ikt(k) e W-{x;k) e w+ (y;k), y > x, 

( d 2 \ _1 f°° 

R w (k)F{x) ={--^ + W{x)-k 2 \ F(x) = J R w (x,(;k)F(0d(. (3.12) 

Remark 3.1. Note that these expressions, originally defined for k ^ 0, are easily 
extended to the point k = 0, for generic potentials. Indeed, one has by Definition al. 0\ 

1 / ; \ , ,\ f-(y,k) f+(x;k) 

e w _{y;k) e w+ {x;k) - 



-2ikt(k) ' w + v ' ' Wron(f + (-,k),f-(-,k)Y 

In the generic case, this expression has a limit when k — > by (|3.9[) and (|3.10l) . 
In the following, we work with the distorted plane waves, which sometimes lead to 
expressions which are only defined for k ^ 0. By the above considerations, it is easy 
to check that in the case of a generic potential, these expressions have a well-defined 
finite limit when k — > 0. 

In particular, we have the following 

PROPOSITION 3.7. Let F e L X (R). Assume W(x) satisfying Hypotheses (V) and 
k G K satisfying Hypothesis (G). Then the inhomogeneous equation 

—^+W(x)-k 2 ) U = F (3.13) 



dx 

has the unique outgoing solution U = Rw(k)F. Moreover, \\U\\l,<x> < C with 
a constant, C(K). 

Proof. Existence follows from the explicit integral representation p,12[) . Note 
that if W is generic, then Rw{k)F is defined for any k £ R, whereas in the non- 
generic case, Wron (/+(x; k), f-(x; k)) — > (k — > 0) and f±(x, k) does not tend to 
zero as k — > [3], so that Rw(k)F has a simple pole at k = 0. 

To prove uniqueness, note that if the difference, d(x), of two solutions is non- 
zero, then d(x) is a non-trivial solution of the scattering resonance problem, that is 
{Hw — k 2 )d = 0, d{x) outgoing at |x| — > oo with scattering resonance energy k 2 £ R. 
However, the scattering resonance energies must satisfy 3(fc 2 ) < 0; see, for exam- 
ple, [TB]. Therefore, d(x) = 0. This completes the proof. □ 
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4. Homogenization / Multiple Scale Perturbation Expansion. 

4.1. Multiple scale expansion. In this section, our goal is to formally obtain 
the expansion displayed in Theorem 12.11 using a systematic two-scale / homogeniza- 
tion perturbation scheme. A proof (and derivation by other means) of this expansion 
is presented in section [5] 

We seek a solution of the equation 

+ y o( x ) + 1 (*> ~) ~ k2 ) e v*+& k) = 0, (4.1) 

in the form of a two-scale function, ey+{x;k) = U e (x, j) which satisfies the jump 
conditions (|3.3[) and the outgoing radiation condition of Definition 13.61 Treating x 
and y as independent variables, we find that U t {x,y) is a solution of 

^ + + V o(x)+q^y)-k 2 ^U*(x,y)=0. (4.2) 

We then formally expand U e {x 1 y) as 

00 

U*(x,y)=J2z j Ufay), (4.3) 
j=o 

and require that 

Uj(x,y + 1) = Uj(x,y), j>0, 

U (x,y) - e lkx , U {x,y) j > 1 outgoing as \x\ -> 00, (4.4) 
Uj(x,y)\ v _ x , satisfies jump conditions (|3.3p . 

The problem is solved by substituting the expansion (|4.3p into (|4.2p and impos- 
ing the equation, jump conditions and radiation condition at each order in e. The 
differential equation becomes 

{-{§- X + \w v ) '' + yo(x)+ q (x,y)-k^j U*(x^)= £ <? r, = 0, (4.5) 
implying the following hierarchy of equations at each order in e 



0(e~ 2 ) 


T-1 = 


-d 2 U 


= 0, 












(4.6a) 


(Die' 1 ) 


r-x = 


-dp! 


- 2d y d x U Q 


= 0, 










(4.6b) 


O(e ) 


ro = 


-dp, 


- 2dyd x U! 


- dpo - 


H(Vb- 


hq)(7 


- k 2 U 


= 0, 


(4.6c) 


Oie 1 ) 


T\ = 


-dp 3 


- 2dyd x U 2 


- dp! - 


~(V - 


Vq)Ui 


- k 2 U! 


= 0, 


(4.6d) 


0(e 2 ) 


T2 = 


-dpi 


- 2dyd X U 3 


- dp2 - 


~(V - 


Vq)U 2 


-k 2 U 2 


= 0, 


(4.6e) 


0(e 3 ) 


»"3 = 


-dp. 


- 2d y d x Ui 


- dp 3 - 


~(v - 


^q)U 3 


-k 2 U 3 


= 0, 


(4.6f) 



0(e*) rj = r, [U j+2 , U j+1 , Uj] = . (4.6g) 

For example, to construct an approximate solution of (|4.2j) satisfying (|4.4[) up to the 
order 3, we solve simultaneously the equations rj — for j = — 2, • • ■ ,3. This will 
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determine the functions Uq, U\, XJ<i and U3 which make U e an approximate solution 
through order 0(e 3 ). Since ey^ + {x; k) — e lkx is to be outgoing, we require U — e lkx 
and each Ui (i = 1, . . . , 3) to satisfy the outgoing condition. We now proceed with 
the implementation. 

Caveat lector! The formal expansion presented in the remainder of this section 
yields terms involving spatial derivatives of ey +(x;fc) and qj(x) of arbitrarily high 
order. Now d x ev +(x;k) has jump discontinuities on supplying and qj(x) has jump 
discontinuities. Hence, the expansion must viewed in a distributional sense, e.g. in- 
volving terms, such as d"d(x — xj) etc. Furthermore, when we impose the jump 
conditions p.3[) to the expansion, order by order in e, we shall throughout assign 
[d x 5{x — Xj)] x=Xj = 0. Although seemingly risky, in section^ we give a complete 
rigorous proof of the expansion with error bounds. 

Beginning at 0(e -2 ), one has from (|4.6a.|) 

r_ 2 = d 2 U = U (x,y) = U (x). (4.7) 
Consequently, one has from (|4.6b[) 

r_i = 3Jl/i = -2d y d x U = Ui(x, y) = U^x). (4.8) 

Recall that y 1— > q(x,y) is 1-periodic and q(x,y) dy = 0. Integration of the 
equation (|4.6cj) with respect to y yields: 

r 1 d 2 

J r (x,y)dy = 0=> --^ U (x) + V (x)U {x) - k 2 U (x) = 0. (4.9) 

Furthermore, since Uq — e %kx is outgoing, one has by Proposition 13.51 

U (x) = e Vo+ (x;k). (4.10) 

By P~9]) and (|4~6c|) leads to 

ro=0 ^ -d 2 U 2 (x,y)+q(x,y)e Vo+ (x) = 0. (4.11) 
Thus, we decompose U2 as: 

U 2 = U^(x) + U^(x,y), 

with U^ p \x,y) a particular solution, and U% (x) an homogeneous solution to be 
determined. 

Again, since y M- q(x,y) is 1-periodic and J Q q(x,y) dy = 0, when by (|4.6dl) . 

f 1 d 2 
J r x {x,y) dy = 0^—^U 1 {x) + V {x)U 1 (x)-k 2 U 1 (x)=0. (4.12) 

Since U\ is outgoing, we claim 

Ui = 0. (4.13) 

Indeed, in this case k 2 is a scattering resonance energy and U± its corresponding mode. 
Scattering resonances necessarily satisfy Qk 2 < [IB]. However, k 2 £ R and hence 
Ux = 0. 
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Consequently, 

ri = =^> -d*U 3 {x,y) - 2d v d x U 2 = 0. (4.14) 

In the same way as for U2, we decompose U 3 as 

U 3 = U^(x) + U^(x,y), 

with U^y (x, y) a particular solution, and C/g (x) an homogeneous solution to be 
determined. 

Integration of the equations (|4.6cp and (|4.6fj) with respect to y, respectively, yields: 

- ^f'W + V {x)ui h \x) + J q(x,y)U^(x,y)dy - k 2 U^ h \x) = 0. (4.15) 

- ^ut\x)+Vo(x)u! i h} (x)~k 2 ut\x) =-J Ui p) (x,y)q(x,y) dy. (4.16) 

We now solve (|4.1ip . (|4.14[) , (|4. 1 5[) and (|4.16|) to obtain a unique (approximate) 
solution satisfying both outgoing and jump conditions, as we see in the following. 
First, we use the decomposition in Fourier series of q(x,y) in y : 

q(x,y) =5>(*)e 2i ™ '. 

Consequently, equation (|4.1ip leads immediately to 

t (p)^ ^ _ e Vo+ (x;k) y, g 3 {x) o2injy 

b'l>i 

From (|4.14|) , one deduces 



U!f\x,y) = - eVo+ ^ 2 ' Y, ^^ 1V - (4.17) 



d 2 U 3 (x,y) = - dx ( ev °+( x; c 2miv_ 
b1>i J 

A particular solution U 3 (x, y) is therefore given by 

Uto\ xiy) = £ ^Me**, (4.18) 

^ b1>i J 

Then, using the Fourier series of q and C/ 2 (p) , we obtain the following equations 
from flUHU) and (|4TT5t : 

" ^^(*) + M^\x) - M h \x) = e -^l £ « and (4.19) 

b'l>i 



-^U^{x) + V,{x)U^\x) - fc 2 [/^(x) 



/>>(*,,),(*,*) * = AE ^ (ey ° +(X;fc) 3 9j(X))9 " 3(X) - (4-20) 
J ° 4?F bl>i J 
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By Proposition 13.71 equations (|4.19[) and (|4.20j) have unique outgoing solutions. We 
refer to the expansion of U e obtained in this way as the 

Bulk (homogenization) expansion: 

U^(x,y) = e Vo+ (x;k) + e 2 (u^ (x,y) + (x)) + e 3 (u { 3 p \x,y) + t/f >(*)) + ... 

(4.21) 

It consists of a leading order average term (homogenization) plus correctors at 
each order in e due to microstructure. 

Failure of Jump conditions at interfaces: 

Recall that we seek a solution which satisfies the jump conditions (|3 . 3[) on U e (x, y) 
for all (x,y) = (x,x/e) at each order in e. The leading order term, ey + satisfies all 



jump conditions. Now consider the terms J7- (x, y) + Uj (x), arising at order Q(e J ). 
By construction, Uj satisfies F 



However U \ (x,x/e) does not. Indeed, for the 
cases j = 2,3, referring to expressions (|4.17[) and f|4. 1 8[) we observe violation of (|3.3I) 
in Uj (x,x/e) at discontinuities of Qj(x) and ey 0+ (x;fc), and their derivatives. 

More precisely, the jump conditions for fail at a; (2 = 1, ... , M) each point 
of discontinuity of q(x, x/e), since one has 



r 2,a' 
e ( - T 2,a 



m 



(4.22) 
(4.23) 



with 



G% 



m 



\3\>l 



— I 
2^ 



2iTrja/e 

3 

e 2i7rja/t 



J 



47T 2 



e Vo+ (a;k) [ qj ] a 

bl>i 

3 i d ^( e V +(x;k)q j (x))} a 



,2iirja/e 



;2 



IjI>i 



(4.24) 



In the same way the jump conditions for fail at points of discontinuity of the 
functions q(x, x/e) and d x q(x, x/e), and for x G {xq, ■ ■ ■ , zjv-i} the support of V S i ng 
( recall: V sing = c i 5{x~Xj) ): 



U<?\x,* e ) 



Tx U ^ 



Ft 



= l -G 



3.u 



H. 



3,aJ 



(4.25) 
(4.26) 



with F 3 a and -ffj a bounded highly oscillating functions and 



G 



3, a. 



2tt 2 



^ l d x( e v +{x\k) qi {x))] a 



2iirja/e 



bl>i 



(4.27) 
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F| a and i/J a can be made explicit, but we omit these expressions as they contribute 
only at 0{e i ). 

Restoring the Jump Conditions at interfaces: 

In order to restore the jump conditions (|3.3[) , we must add to the expansion, 
at each point where the jump conditions are not satisfied, an appropriate corrector. 
These correctors each solve a non-homogeneous equation, driven by the jumps in the 
bulk expansion (|4.21[) . 



To see this, first note that 



d_jAp) 

dx U j 



©(e" 1 ), j = 2,3. Since u\ p) contributes 



at order e J , this suggests adding a corrector at order e 3 . Thus, we introduce the 
Bulk expansion with corrector terms: 

U\x,y) = e Vo+ (x;k) + eU{{x) + e 2 ( U%>\x, y) + C/f \x) + «!(*)) 

+ e 3 ( U^\x,y) + Ut\x)+Ul(x)) + ... (4.28) 

The interface correctors Uj(x) are to be determined so that, at each order in e, the 
expansion (|4.28l) satisfies the jump conditions (|3.3|) . the differential equation (14.11) 
and outgoing radiation condition. 

We construct Uj(x), j = 1,2 below. The general construction uses the following 

Lemma 4.1. Let Fi,F 2 G E and V = V smg + V reg as in (|2.1I) . Then there exists 
U(x), an outgoing piecewise C 2 solution of: 

d 2 \ 
-— + V (x) - k 2 \ U = 0, for x < a and x > a, (4.29) 

which also satisfies the following jump conditions at the point x = a: 

P{x)] a = F U 

- cU{a-) = F a . 



ax 



Here, U(a—) = VnnU(x), and the constant 



if a £ supp V sing , 
if a = x jo e supp Vsmg ; 



(4.30) 



recall V smg (x) = Y^f^ c 3 S(x-xj). 
U(x) has the form 



u{x)={ rM* ; S 'i x<a ^ (4.3i) 

v ' \ f3e Vo+ (x; k) if x > a, y ' 

for appropriate choice of a and j3, namely 

F 2 e Vo+ (a;k) - F 1 d x e Vo+ (a+;k) F 2 e Vo ^(a;k) - F 1 d x e Vo ^(a+;k) 

a = ; — and a = r . 

2ik t^ orn {k) 2ik t^ om (k) 

(4.32) 
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Before giving the proof, we explain why choosing Uj as in Lemma 14.11 does not 
change the bulk expansion (|4.21[) constructed above. Therefore, our approach which 
first computes the bulk-expansion and then the correctors is consistent. 

As pointed out the expressions in the bulk expansion (|4.21[) 



U {x,y) = U^ lk (x,y) = Uf> (x, x/e) + U\ n > (x) 



(p), 



(4.33) 



do not satisfy jump conditions (|3.3p . Suppose now that we replace the functions 
Uj(x,y) = U!j ulk {x,y) by U,{x,y) = U^ lk (x, y) + U e a {x) and we seek U e a (x) so as to 
ensure jump conditions (|3 . 3[) . (Assume only one corrector is required). Note that 
since Uq{x) lies in the kernel of d v , adding such a term has no effect on the equations 

determining Uj P \x,y). Further, we want to preserve the form of Uj h \x), which has 
previously been constructed. Thus, 

r 3 [U j+2 , U j+1 , Vf lk + W a ] = rj [U j+2 , U j+1 , U^ k ] 

+ (-d 2 x + V {x)-k 2 )W a {x) + q(x,y)K(x). (4.34) 

The equation for Uj h \x) is obtained by averaging (|4.34[) with respect to y. Since 
q(x, y) has mean zero with respect to y, this gives 



rj [U j+ 2,U j+u Uj 



bulkl 



(x,y)dy + (-d 2 x + V {x)-k 2 )Ul(x) = 0. (4.35) 



Thus, if we choose U^{x) to satisfy (|4.29[) . then the second term in f|4. 35[) vanishes and 
the equation for Uj (x) is preserved. Therefore, if Lemma 14. II is used to determine 
the jump-driven correctors at each order in e, then the corrected bulk expansion (|4.28|) 
is the solution we seek. 



Proof of Lemma \4-l\ The piecewise form of U (|4.31|) satisfies the outgoing radiation 
condition, by construction. The constants a and ft are determined by the jump 
conditions. 

Using the fact that ev +(x;k) and ev -{x;k) satisfy the jump conditions (|3.3[) , 
one has 



[U{x)] a 



[±U{x)} a - oU(c 



(3e Vo+ (a; k) - aey _(a; k), 

j3d x e Vo +{a+] k) - ad x e Vo -{a-; k) - ca e Vo -(a—; k) 
(3d x ev +(a+; k) - ad x ev -(a+; k). 



Solving this inhomogeneous system, using the value of the Wronskian, given in (|3.8I) . 
leads immediately to (|4.32[) . This completes the proof of Lemma 14.11 □ 

We now proceed to apply Lemma 14. II to determine the correctors associated with 
U^ p) and Using (|4~22]) - (|4~23|) and (|L25)) - (j4T26)) . the jump conditions (O applied 
to eU{ 



2 TT {P) 



e z U. 



e 3 U^ p) read: 



dx Ul 



e [Wf]„ + c 2 ( Fl a + [W|]„) = G(e 3 ), 
a -dA{{a-)) + e 2 (%- C f/f(a-) + GL 

ax 



-cU{(a-))=0{e s ) 



(4.36) 



(4.37) 
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Equations (|4.36|) and (|4.37|) imply jump conditions at order e and order e 2 . There- 
fore, we construct Uj a , j — 1,2 solving the two inhomogeneous problems at each 
point, a, of non-smoothness. 
System for corrector U{ a : 



(4.38) 
(4.39) 

(4.40) 



+ y ( x )_fc2j W e o=0i x ^ a , 



[Ka] a = 0, 



dx 



-G% 



System for corrector W| a : 



-F 



2,ai 



dx 



Ula{x) -cUl a {a-) = -Hl a -Gl a +cU^\a-). 

(4.41) 

Lemma l4~Tl applied to (|4.38|) - (|4.39[) and (|4.40p - (|4.41[) defines the unique correctors 
U{ a and W| )0 : U\ a is given by (j4.31[) . i.e. 



Ut a (x) 



a\ a ev ~(x;k) if x < a, 
Pt,a e v +(x;k) if x > a, 



(4.42) 



with a\ a and {5\ a given by 

G%a 



~2ikt% om (k) 



e Vo+ (a;k), Pl <a = 



G 



— — e Vn ^(a:k), (4.43) 
2ikt>> om (k) y ' ; ' V ; 

where G\ a is given in (|4.24j) . Then, a is given by (|4.31|) with a\ a and /3| a a and 
ft given by 



^ om (fc) (( 



Via = 7T- 



2ikt% om (k) 
1 



2ikt% om (k) 



Ka ~ Gla + c U^(a-)) e Vo+ {a- k) + Fl a d x e Vo+ (a+; k)) , 

a + c ^ p) (a-)) ey _(a; fc) + *£ o x e Vo _(a+; fc)) , 

(4.44) 



where fff a , F| a and G 3 a are given in (|424)l and (|4~27)) . 
Therefore at 0(e), we define the corrector U{ as 



M 



(4.45) 



where Oj, j = 1, . . . , M denote the points of discontinuity of q (x, x/e) . 
At order C(e 2 ), we have a violation of the jump conditions (|3.3j) due to 

(i) points of "discontinuity" of qj(x), i.e. aj (j = 1, . . . , M) for which [qj] a ^ 
or [d x qj] a ^ 0, and 

(ii) the singular set supp V S i ng = {xq, ■ ■ ■ , xn-i}- 
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Thus we construct, U\ a , for all a in the set, f2, of non-smooth points of V e (x): 
Q = {x , x N -i} U {-00 = a , ax, . . . , a M = 00} (4.46) 
and define the corrector t/J by: 

Z4 = E^V (4-47) 

We summarize the preceding calculation in the following 
Proposition 4.2. 

ey+(x;k) = U e (x, x/e) = e Vo +{x; k) + eU{(x) 

+ e 2 ( U^ h \x) + U<*\x,x/e) + W 2 {x) ) +0(e 3 ) (4.48) 

gives a formal construction of the distorted plane wave ey + {x;k), through 0(e 2 ) 
with error of size C(e 3 ). The correctors U{(x) andti^ix) are given by (|4.38[) - (|4.39f) 
and (|4T40)) - (jOT|) . 

Finally, u!j p \x,y) andU^ h \x) are given by 



4?r 2 ^ «v~/ -2 



(-£ + %(*) - fc 2 ) E^(*) = £ ^ ^ Agoing, 

W'(x,y) = -—3 2^ 9x(evb+(a:;*!)?j(x))— ^3— , 
*" IjI>i 3 

~ + Vo{z) - kA Ui h \x) = i J2 d -^Y^0MhzM , uW outgoing . 
X ' b'l>i n 3 

4.2. Expansion of the transmission coefficient, t £ (k). The results of the 
previous section can now be used to derive expansion ()2.13|) for the transmission coef- 
ficient, t e (k) associated with the potential V e (x). t e (k), through order e 2 is derived by 
isolating appropriate terms in the expansion (|4.48[) . The sense in which the remainder 
is small is proved, by entirely different means, in section [5] 

0(e ): The only term at order one is ey +, which gives the leading order transmission 
coefficient, i lom (fc), corresponding to the average potential Vo- 

0(e 1 ): At order e, we seek the contribution to t t (k) from U\. From (|4.31l) we have, 
since ey 0+ (x; k) ~ ig om (fc) e lkx as x — > +00, that the contribution of U{ a to 
the transmission coefficient is given by 

1 „2iirja/e 

*!,«= Pt, a to° m (k) = ^e Vo+ (a;k)e Vo -(a;k)Y, [<&]„ ~~ ■ 

\i\>i 3 

Finally, summing over all contributions from points of discontinuity of qj , one 
obtains the complete first order contribution from IA{: 



ST e Vo+(Q3; fc )evo-(Qj; fc ) , , 
h - 2^*i.«i - 2- 4^1 1^ — J — ■ ^ A9 > 

3 3=1 \1\>1 
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C(e 2 ): (a) No contribution to i e (fc) from e 2 Jjif^ : We estimate pointwise 



1^ (*, f) I < ^|evb+(*;*)l E T ^ C ( E 



bl>i 3 b'l>i 



-s> 0, x oo. 



Here we have used the uniform bound (|B.3p on ey 0+ for x > and the 
hypothesis (|2.9p . Since, fTg (x, | ) — >■ as x — > oo, it does not contribute to 
the transmission coefficient. 

(b) Contribution ofe 2 U^ h) (x) tot e (k): 
From (|4.19[) . one has 

7 |j|>i J 

Using expression (|3.12p for the outgoing resolvent we have 

\ lil>i 
-1 [ x e Vo+ (C;k) ^ \ qj (0\ 



f x ey 0+ (C;fc) v \ qj (0\ 2 ( , , v , , v 



'2ikt h ° m {k) J_ x 47T 



2ifc tQ° m (fc) Jj. 



bi>i J 

Therefore, since qj £ L 2 for all j £ Z, and ey ± € L°°, one has when x — > oo: 
/ -ey +(a;;fc) /" ^ \qj(0? 



It follows that the contribution of U^ix) to the transmission coefficient is 

'(*) = --^ / E ^^W -(C; fc)ey 0+ (C; *K- (4.50) 

J R J 



±hom ( 
l 2 



bl>i 

(c) Contribution of e 2 £/J t e (k): 

We study as above. From (|4.31[) we have, since ey 0+ (x; fc) ~ t^ om (k) e lkx 
as x — > +oo, that the contribution of £/J to the transmission coefficient is 
given by t% = t%° m {k) 
From (|Q3|) . P~2"4"|) and (|43d| we have 

*!,»(*) = 0!,a 

= J g^" 2 . fc — E ( c e vb+( a ; k)qj(a-) + [d x (e Vo+ (x; k) qj {x))] a ) - — 

n % b'l>i J 



irr E e Vb+(a; fc ) fe] a ^ey _(a+;fc) 



&7T 2 ik ^-"tn-n^— -v-rv . 2 
lil>l J 
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Using the easily verified identity 

[ qj(x)d x e Vo -(x;k) ] a = d x e Vo -(a+;k) [ qj(x) } a + c e Vo -(a;k) qj(a-), 

(4.51) 

we obtain 

^'"^ = s^ikYl ([4(eF»+(i;%W)] a evo-(«;fc) 

bl>i 

e 2inja/e 



[qj(x)d x e Vo -(x; k)] a e Vo+ (a; k) 



J 2 



, }2 l d ^( e Vo~(x;k)e Vo+ (x- : k)q j (x))] a — ^ — 



b'l>i J 

Finally summing over all contributions of all singular and / or discontinuity 
points of V e , we obtain the simple expression: 



t 



2%-Klaje 



2( k ) = XX«0) = -^y, i d Aev ~{x;k)e Vo+ {x;k)qi{x))] a — - — . 

aeQ agn|Z|>l 

(4.52) 

C(e 3 ) : By similar considerations to the above discussion of and , the terms 
[/J = [7g + in the expansion of ey^ give a correction to t e {k) of order 
e 3 , and is therefore subsumed by the error term in the expansion (|2.13|) . 

In summary, we have an expansion of t e (k), agreeing with the expansion (|2.13l) 
in Theorem 12.11 

Proposition 4.3. Formal corrected homogenization expansion: 

t e (k)=t$ om (k) + et\(k) + e 2 (t'i om (k)+t' 2 (k)) + 0(e 3 ), (4.53) 

where the leading order term, tQ OTn (k), is the transmission coefficient associated with 
the homogenized (average with respect to the fast scale) potential Vq, t2° m (k) is a 
classical homogenization theory corrector given by (|4.50|) . andtj, j = 1,2 are interface 
correctors given by (|4.49p and (|4.52|) . 

Note that if Vq is generic, then since using that ty (k) and ev ±(x,k) are O(k) 
as k — > 0, we see the expansion is formally valid for any However, ifVo is not 

generic then we must exclude k = 0; see the discussion of Remark ] 3. 11 

5. Rigorous analysis of the scattering problem. In the preceding section, 
we applied the classical method of multiple scales to derive a formal expansion for the 
distorted plane wave ey+(x; k) and transmission coefficient t e (k); see section[3] For 
sufficiently smooth potentials, this expansion satisfies, at each order in e, all necessary 
continuity conditions as well as the radiation condition at infinity; see Definitions 13. II 
and EH 

We found, however, that if the potential is non-smooth this expansion, while valid 
in the bulk, violates continuity conditions at (i) discontinuities, and (ii) at strong 
singularities of the background, unperturbed potential, Vq(x) = V reg (x) + V S i ng (x). 
We found, in Section!?] that we can, "by hand", construct interface correctors for each 
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point of non-smoothness, thereby giving a corrected expansion (bulk expansion plus 
interface correctors) which is a valid solution to any finite order in e. The expansion 
of Proposition ^. 31 is explicit through order e 2 with order e 3 correctors. 
Question: Does the procedure of section [4] yield a valid expansion with an error 
terms satisfying an appropriate higher order error bound? 

It turns out that the formal expansion is correct with an appropriate error esti- 
mate. However, we obtain this result, not by expansion in scalar e but rather in the 
the function q e (x), with respect to which there is an analytic perturbation theory in 
an appropriate function space . Smallness required for control of the perturbation ex- 
pansion derives from q e (x, x/e) being supported at high frequencies if e is small. The 
principle terms, displayed in the expansion of Proposition 14.31 (and indeed the terms 
at any finite order in the small parameter, e), are obtained via small e asymptotics of 
the leading order terms in the q e expansion. The approach we use was introduced by 
Golowich and Weinstcin in [9]. 

5.1. Formulation of the problem. We consider the general one-dimensional 
scattering problem 

(—^2+ Vo(x) + Q(x) - k 2 ^ e v+ = 

e Vo+Qt+ (x;k) - e ikx -> 0, x -> -oo, (5.1) 

where Vo (x) as hypothesized in section [T] and Q is a spatially localized perturbing 
potential, which we think of as being spectrally supported at high frequencies. Q may 
be large in L°°. As a model, we have in mind Q(x) = q e (x) = q{x,x/e), with e small. 
We introduce the scattered field, u s , via 

e v +(x; k) = e Vo +(x; k) + u s (x; k) (5.2) 

where u s is outgoing as x — > ±oo. Therefore, u s is the solution of 

+ ( y o + Q)( x ) - fc2 ) U ^ k ) = -Qzv a +{x\ k), (5.3) 

with outgoing conditions : ( d x =p ik ) u s — > 0, x — > ±oo. 

Applying the outgoing resolvent, Ry (k), to (|5.3[) and rearranging terms we obtain 
the Lippman-Schwingcr equation 

u s = - (I + i?y (fc)Q) _1 R Vo (k) Q e Vo+ (-;k) =>> 
e v+ (x;k)= e Vo+ (x;k) - (I + R Vo (fc)Q)" 1 R Vo (k) Q e Va+ (-;k). (5.4) 

Consider now the formal Neumann expansion, obtained from (|5 .4[) . 

e v+ (x; k) = e Vo+ (x;k) - R Va (k)(Qe Va+ (x; fc)) + R Vo (k)QR Vo (k)(Qe Vo+ (x; k)) + . . . 

OO 

= e Vo+ (x;k) + ^2 l(- R v (k) Q) m e Vo+ ](x;k). (5.5) 

m=0 

In this section we show for a class of Q, which include high-contrast (pointwise large) 
microstructure (highly oscillatory) potentials that the expansion (|5.5[) converges in an 
appropriate sense and that any truncation satisfies an error bound. 
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5.2. Reformulation of the Lippman-Schwinger equation and the norm 

IIIQUI . We seek a reformulation of the Lippman-Schwinger equation (|5.4[) in which it 
is explicitly clear that if Q is highly oscillatory, then the terms of the Neumann series 
are successively smaller. Introduce, via the Fourier transform, the operator (Dq) s 

i r +oa 

(D ) s g= (I- Ay/ 2 g ^—J e«« (1 + £ 2 ) s/2 (5.6) 

and the localized function X 

X (x) = (x)-° = (l + x 2 )" f , a>4. (5.7) 

Now introduce the spatially and frequency weighted distorted plane wave, Ey+(x] k), 
given by: 

E v+ {x-k) = ( (D ) X e v+ )(x;k). (5.8) 

With the operator definitions 

T Rvg (k) = (D ) xRv (k) X (Do) , (5.9) 
T Q = (Do)- 1 X^QX- 1 (Do)' 1 = (Do)' 1 (x)° ■ Q ■ (x)° (Do)' 1 , (5.10) 

Ey+(x] k) can be seen to satisfy 

(/ + T RVq Tq) ( Ev+{-\ k) - E Vo+ {-- k)) = - (Do) xRvM Q e Vo+ (x; k). (5.11) 

Here's the motivation for our strategy. Note that Tq has the operator (Do) 1 
as both a pre- and post- multiplier. This has the effect of a high frequency cutoff. 
Therefore, for highly oscillatory Q, Tq is expected to be of small operator norm. If 
the norm of Tr v o Tq is small then / + Tr v o Tq is invertible and we have the 
preconditioned Lippman-Schwinger equation 

E v+ (-;k) = E Vo+ (-;k) - (7 + T^Tq)" 1 (D ) X Rv {k) Qe Vo+ (x;k). (5.12) 

We proceed now to construct a norm, |||Q|||, such that if |||Q||| is small then 
Tr v Tq is bounded and of small norm as an operator norm from L 2 to L 2 . 
The norm we choose for the perturbing potential is defined as follows: 



\Tq\\ L 2_> L 2 



[DoT 1 (xYQ (xY (Do)' 1 l2 ^ l2 , a > 4. (5.13) 



The next result establishes the expansion of the distorted plane waves ey+(x\ k) in 
a iJ 1 (IR; x{x)dx) and therefore, by the Sobolev inequality, a L°°(R; x(x)dx) convergent 
expansion for |||Q||| sufficiently small. 



Theorem 5.1. Let V satisfy Hypotheses (V), and k € K a compact subset of I 
satisfying Hypothesis (G). Define 

t (K) = n 1 - -I, > 0. 

maxfcgK \\TR VQ (k)\\ L2 ^ L2 

I f\\\Q\\\< to{K), then for all k £ K: 
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• The preconditioned Lippman Schwinger equation (|5.12[) has a unique spatially 
and spectrally weighted distorted plane solution, Ey + (x; k). 

• This solution can be expressed as a series, which converges in L 2 (K), uni- 
formly in k G K: 

oo 

E v+ {x-k) = E Vo+ (x;k) + J2 {-T Rvo (k)T Q ) m [(D )xRv (k) Q e Vo+ ](x 

m=0 

= E Vo+ (x;k) - (D Q ) x Rv (k) Q e Vo +(x;k) 

+ T Rv ( k ) T Q ( D o) X Rv a (k) Q e Vo+ (x; k) - . . . 

• It follows that the distorted plane wave, ev +Q,+ (x;k) satisfies the approxi- 
mation for any M > 1 



[Do)x (u.(-,k) + J2 [(-Rv (k)Q) m e Vo+ ](-,k) J 

\ rn=0 / 



< (5.14) 



withu s (x;k) = e Vo+Qt+ (x;k) - e Vo+ (x;k). 



Remark 5.1. In the proof of Theorem I5.il the distinction between generic and 
non-generic cases arises through the properties of the unperturbed resolvent, Rv (k), 
as k — > 0; see Proposition \3.7\ 

In the following, we prove that both Tr v and Tq are well-defined operators, 
bounded in L 2 . Then, Theorem 15.11 follows immediately if Q satisfies the smallness 
condition 

\\ t q\\l^ L 2 < j™g (ll T ^ n ( fc )|| L 2^ L 2) _1 = r (k). (5.15) 

Proposition 5.2. Let (x) 2a Q(x) e i 2 (R). Then T Q , as defined in (|5.10j) . is a 
Hilbert- Schmidt operator and is therefore compact. 

Proposition 5.3. Let q e (x) satisfy the conditions in Hypotheses (V). Then, for 
e small, 



Proposition 5.4. T Rva (k) is a bounded operator from L 2 to L 2 . 
Propositions 15.21 and 15.31 are proved below. The proof of Proposition 15 .41 is some- 
what more technical proof and is found in Appendix [Cj 

We now prove Proposition [521 The proof of Proposition 15.31 begins on page [231 
Proof. Proof of Proposition 15. 2t We begin by introducing the notation 



Q^x-'Qx- 1 - 



(5.17) 
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Then, one uses the following calculation: 
(A))" 1 Q« (Do)- 1 f(x) = (Do)- 1 Qt (Do)- 1 (J- J e **/(0 d£ 

c 47T 



^3 (| e-^/(C) dCj J mWr^&it-v) dr, d£ 



f(OK(x,C) d(. 



C 



with the kernel 



k( X , o = ^JJ(i+ er 1/2 (i + ify 1,2 z i(xv - a) QKt - n) d V (s.is) 

We want to prove that fj \K(x, C)| 2 dx d( < +00, i.e. K e L 2 (M 2 ). One has 



K(s, z)= 1 1 K{x, C,)e~ lxs e~ Kz dx dQ 

(0' 1 (v)~ 1 e ix(r '~ s) e- i ^ + ^Qi(^-'n) dx d( drj d£ 



1 

1 QH-s-z) 

~87T 3 (l + S 2 ) 1 /2(l + z 2)l/2' 

Therefore, we deduce 



Since Q" e £ 2 (R), one has immediately ff xl - \K(x, C)\ 2 dx d( < 00, and 

H-^C^j 1| z, 2 (M 2 ) — ^ll^ lll,2(K)- 
Therefore Tq is a Hilbert-Schmidt integral operator, and is therefore bounded, with 

II -^Q II i=— >.£= — ^1|Q"IIi,2(r)- 

This completes the proof of Proposition 15.21 □ 
Proof. Proof of Proposition 15. 3t 

Consider Tq, where Q = q e (x) — q (x, -) as in (|2.2I) . From the proof below, one 
has T qc f(x) = f(£)K e (x,£) d(, with the kernel K e (x,() satisfying 

\ 1 ql(-s-z) 
K e (s,z) = 



87T 3 (1 + S 2 )V2(1 + z 2)1/2' 
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Using the decomposition in Fourier series of q e (x) = q (x, -) , one has 

q\{x) = X^qeX- 1 ^) = E 4(^ 3{x/e \ 
lil>i 

and therefore 

S(o = E / ^y-^e-^dx = e 5 (— - 1 
\j\>i Jx u\>i v e 

One deduces then 



|,,!( '' +2 ' ^3/e),l 

, V- /" , |^ + (2^/ e ))| 2 



S^^Vk 7 r '(1 + (7 ? -z) 2 )(1 + z 2 ) 

|g(r ? + (27rj7e))| 2 
H >^ "(l + (»7-z) 2 )(l + z 2 ) 



Estimation of Ji: 



l?5(r? + ^)| 2 



drj 



h.A = 



W >st (l + ( V -z) 2 )(l + z 2 ) 
h,A + h,B, with 

dz [ di]<C% 



dz 



dz 



M> 2 



$ h + — ) 



| 2 |>ii l + fa-z) 



drj 



;dzdr] 



< C ^\\q)\\U ^d as \ V - z\ > ^ for \z\ < \ V \ > ^, 



"J 



2c 



2r 



h,B = 



dz 



drj<C-^ 



\n\> 2 



F J\»\<% 1 + ^ 



dz 



c 



drj 



< & 



Now, summing on j, one obtains I\{e) = 0(e 2 ). 

Estimation of I^'. We first show that if we assume only that 5Zh'|>i Ikjlli 2 < °°' 
then h{e) -»• as e ->• 0, and therefore ||T g J|| 2 ^ L2 = o(l) + C(e 2 ) = o(l) as e -> 
with no specified rate. 

We then show that if q t is as in hypotheses (V) then \\T qc \\l 2 ^,l 2 — C(e) as e — >• 0. 
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Assume J2\j\>i hjWh < oo. Then, 



h = E / dz 



Igjjfr +(27Tj/ e ))| 2 

H<^ ?? (l + (r ? -^) 2 )(l + ^) 



< 



(l + z 2 ) 



bl>i 



|d(r ? +(2 7 ri/e))| 2 dr 7 



* C E 

b1>i 



l<7*(r)| 2 dr. 



Note that £ m>1 Jr kS( T )| 2 ^ r = Eh|>i WljWh < °o, implying I 2 = o(l) as e -> 0. 



We now turn to the case where q e satisfies the condition in Hypotheses (V) in 
order to establish that liT^JU 2 ^ 2 = 0(e) as e — > 0. The estimate for 7i(e) is as 
above: h{e) = 0(e 2 ). 

Now, we estimate /2(e) using the fact that since g!j <E L 2 and (q^)' € L 2 : 



M 



a.1 + 1 



E / Qj(^- tTX dx 
1=0 Ja i 



M 

1 E 



1=0 KJai 



at + i 



(qh'(x)e- iTX dx - g»(ar +1 )e- ?Tai+1 + d( a '+K 



< C- 









-a) 


a j ) 





Therefore, one has 



h = E / dz 



| g »(r /+ (27rj7e))| 2 
|„|<si , '(l + (l^) 2 )(l + ^) 



< 



One deduces finally that ||T ge || x 2 — >z, 2 = 11^11^2^2) = h +^2 = This completes 

the proof of Proposition 15. 31 □ 



5.3. Application to the transmission coefficient, t(k) = t[k;Q}. This sec- 
tion is devoted to the proof of Theorem 12.21 The heart of the matter is to view t(k) 
as a functional of the perturbing microstructure potential, Q(x) 



Q ^ t[Q] 



(5.19) 



and to use the Lippman-Schwingcr expansion of Theorcm l5.1l to expand t[Q] for small 
IIIQIII: 



t[Q] = C m + h[Q] + t 2 [Q,Q] + t 3 [Q,Q,Q] 



(5.20) 
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where tj [Q, Q, . . . , Q] is j— linear in Q. The transmission coefficient expansion of 
Theorem 12.21 is recovered from the small |||Q||| asymptotics of the first several terms 
of the expansion of t[q e ]. Finally, the error terms are estimated. 

Recall that from (|3.7p the transmission coefficient, t\y(k), associated with the 
distorted plane wave ew+(%\ k), is given by 



We denote the transmission coefficients of ev +(x; k) and ey 0+ Q^ + (x\ fc), respectively, 



tv (fc) = t (k) = C m (k), 

t v (k) ee t(k) = lim e~ ikx evb+Q,+(a:;A) = t'> om (k) + lim e~ ikx u s (x; k). 



To obtain the desired leading order expansion of t{k) of Theorem 12.21 we now derive 
the small |||Q||| asymptotics of the linear and quadratic terms in Q of (|5. 14|) . 

Calculation of ti[Q]: 

One has from (|3 . that 



tw(k) = lim e 



ikx 



e w +(x;k). 




h[Q] e 



ikx 



as x — )■ oo , 



where 




(5.21) 



Calculation of t2[Q, Q]: 



One has from (|3 . llj) that 



Rv {k) Q Rvo(k) Q(0 e Vo+ (x;k) 




HQ, Q] e 



ikx 



where 




(C;fc))ey _(C;fc) d( 




Q(C)ey _(C; k) ( 1,(0 + 7 P (0 ) d(, with (5.22) 



HO = J-ooQ( z ) e v +(z;k) e Vo -(z;k) e Vo +((;k) dz, 
7 r(C) = Q(z) e Va+ {z]k) e Vo+ (z;k) e Va _{C;k) dz. 
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Estimation of the error terms: 

The final step for the proof of Theorem [2]2] consists in a bound on the contribution 
to the transmission coefficient from the remainder term in expansion (|5.20p . This is 
given by the following Theorem: 

Theorem 5.5. Let K denote a compact subset o/M, satisfying Hypothesis (G). 
Introduce for k £ K 

t rem (k;Q) = t(k;Q) - t h a ° m (k) - t x [Q] - t 2 [Q,Q}. (5.23) 

Then we have, uniformly in k € K: 

1. If V has compact support, then t rem (k) = 0(\\\Q\\\ 3 ). 

2. IfV is exponentially decreasing, then t re m(k) = 0(|||<3||| )■ 

3. If (x) p+1 V (x) € L X (R) and (x) p Q(x) <E L 2 (R), p > 8, then there exists 
2 < (3 < 3 such that t rem {k) = O(\\\Q\\\ ). 

Proof. It is convenient to first introduce 

f rem = - (i + Tr^Tq)' 1 (T Rv T Q f (D Q ) xey 0+ (z;fc) (5.24) 
= (Dq)xu s + (D )xRv a (k) Q(x) e Vo+ (x;k) 

- (D ) X Rv (k) Q RvM Q{x) e Vo+ {x;k). (5.25) 

Using (|B.4[) . one deduces that (Do) x e v a +( x 'i G L 2 , with 

||(D )xe\/ +|| L 2 = ||(?7)xe^:(77;/c)|| 2 

V 

<\\x{x)e Vo + (x;k)\\ L2 + \\d x {x(x)e Vo+ {x;k))\\ L2 < \\ X {x)(x)\\ L2 . 

X XX 

Therefore, thanks to Propositions 15.21 and 15. 4[ and using (|5.24j) . one has for |||Q||| 
small enough, 

/™£i 2 and ||/ rem || L2 < C x IIIQIH 3 . (5.26) 
The following pointwise bound can be also be deduced 

(Do) 1 frer 

which implies 

X -\x)(D Q )- 1 f rem {x)\<C x x-'ix) IIIQIH 3 . (5.27) 

From (|5.25[) we have that t rem {k) is the complex number for which 

lim ( X ~Hx) (Dor 1 frem ~ t rem (k)e' kx ) = 0. 

We now use the decay properties of the potential V to estimate the magnitude of 
t r em(k) for |||Q||| small. 

Case 1: V has compact support 

Assume suppU C [—M,M], M > 0. Using the explicit representation of Ry , (|3.12[) . 



< 



(ri)~ l frem{r{)^ X 



< ll( ? /> 1 \\hl ll/rem|| L 2 , 
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for x > M we have: 

1 f x 

-R Vo (k)Q e Vo+ (x\k) = — / Q(Q e Vo+ ((; k) e Vo - (C; k)e Vo+ (x; k) d( 



Q(C) ev +(C;fc) e Vo+ (C;k)e Vo _(x;k) d( 



1 



2z \z tQ 



lio-ni 



+ 00 

<2(C) ey 0+ (C;fc) ev _(C;A;)ey 0+ (x;fc) tiC 



|g}e Vo+ 0r;fc) = h[Q]e ikx 



'() 

Similarly, for the quadratic in Q-term we have 

Rv (k) Q RvM Q e Vo+ (x;k) = t 2 [Q,Q] e lkx . 

Therefore, 

u s = h[Q] e ikx + t 2 [Q,Q] e ikx + t rem e lkx , 

where for x > M 

X" 1 (-Do)" 1 frem{x) = t rem {k)e ikx . 
Therefore, using the pointwise bound (|5.27[) we have 

|W*)I < C xX ~\M) \\\Q\f = 0(\\\Q\\n 

Case 2: V is exponentially decreasing 

Assume |V (a:)| + \Q(x)\ < Ce~ a ^ for some C, a > and x > M. As in the first 
case, the formula for the resolvent (|3.12p leads to 

-R Vo (k)Qe Vo+ (x;k) = h[Q] e lkx + *H (e Vo+ (x; k) - t h Q ° m e ihx ) 



1 



+00 

<2(C)ey +(C; k ) (ev +(C; k)e Vo -(x] k) - e Vo+ (x;k)e Vo -((;k)) d(. 



Onh fhom 

611% Lq j x 

Using (|B.4[) . one can easily bound for x > M 
1 f + °° 

2ik t h om J Q(Q e Vo+(C,k) {e Vo+ {C,k)e Vo -{x;k)-e Vo+ {x;k)e Vo ^((;k)) d( 

r+oo 

<C Ce~ a ^(() d( < C'e~ a/2lxl . 

J X 

Now, we use the estimate (|B.3|) 

\m+(x; k)-l\< 1+ J (1 + \s\)\V (s)\ds, 

so that |ey 0+ (x; k) - t>™ m e lkx \ < Ce~ Q / 2x for x > M. Finally, one obtains 
\R Vo {k)Q e Vo+ (x;k) -t 1 [Q]e lkx \ < Ce- a ' 2x . 
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A similar estimate holds for u s = ey + (x; k) — ev +(x; k), and for the quadratic 
term Ry a {k) Q Ry Q (k) Q ev +{x; k). Therefore, for x > M we have 



X 1 (Do) 1 f remix) - t rem (k)e 



ikx 



< Ce~ a/2x . 



Again the pointwise bound (|5.27[) implies, for x > M, that 



|frem(fc)| < C x \~\i 



+ Ce~ a/2x . 



Finally, choosing x — — — In |||Q|||, one has for |||Q||| small enough, 

IWfc)l< qilQ||| 3 (i + (in|||Q|||»<q||Q||| 3 -. 

Case 3: (x)P +1 V (x) G L X (R) and (x)PQ(x) G L 2 (R), with p > 8 
We use again the formula of the resolvent (|3.12p : 

-R Vo (k)Qe Vo+ (x;k) = t x [Q] e lkx + *-M (e Vo+ (x; k) - t%° m e ikx ) 



1 



2ik t% om 



+ 00 



Q(Oev +((;k)(e Vo+ ((;k)e Vo -{x;k) - ev +(x;fc)ev -(C;fe)J d(. 



Using the estimate (|B.3|) leads to 
\e Vo+ (x;k)-t%° m (k)e ikx \<C 
Therefore, one has 



^-L^(i+\ s \y+i\v (s)\ ds < c^\\v \\ L ,, P+1 . 



i 



2ik t% om 



+00 



Q(0 ev +(C;fc) (ev +(C;k)ev -(x;k) ~ e Vo +{x;k)e Vo -(C;k)) d( 



< C 

from which we deduce 



\RvMQ ev +(x;k) -h[Q] e* kx < 



J0_ 

(x)P 



C 

W 



d( < 



C 

w 



1 



( C ),-1 



10(011 



I Vol 



(cy- 1 



\Q\ 



L' 2 



L- 



Similar estimates hold for u s = ey+(x] k) — ev +(x] k), and for the quadratic term 
Rv (k) Q Rv (k) Q ev +{x; k). Therefore, for x > M, one has 

C 



X 1 (D ) 1 f remix) ~ t rem (k)e 



Ikx 



< 



(X)P 



I Vol 



L 1 :P+ 1 - 



Since x( x ) = ( x ) a with a > 4, the pointwise bound (|5.27[) yields 

1 



trem(k) 5; C*x(^) 



C- 



so that choosing x = \\\Q\\\ 3 ^ p+a \ which tends to infinity as |||Q||| tends to 0, one 
has 



3p 
p + a 



\trem(k)\ < C 

It follows that with a > 4 and p > 2a, one has 

\trem(k)\ = 0(|||Q|||^), 2 < /3 

This completes the proof. □ 



_3p_ 

p + a 
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5.4. Completion of the proof of Theorem 12.11 In this section, we show 
how to derive the corrected multi-scale / homogenization expansion of section from 
the rigorous results of the previous section with a potential V = Vq + q c satisfying 
Hypotheses (V) , and using Proposition 15.31 Theorem 12.11 follows then as a direct 
consequence. 

The small e asymptotics of ti[<7 £ ]: 

We use the decomposition of q e in Fourier series in y 

q e {x)=q[x,-) = ]T qj {x)e 2i ^ x '*\ 

e m>i 

that we plug into t\ [q e ] , given in (15.211) : 

1 f +oc 
kfe] = 2ikT< with = / <z,(C)ev 0+ (C;fc)ey _(C;fc)e 2 ^ (c/e) d(. 

1 b1>i 

We assume that qj is piecewise C 3 , so that there exists — oo = do < ai < • • • < aju < 
oa/+i = oo, such that qj £ C 3 (ai,ai+i). Then, one has 

4'%] = ^^ ai+1 <7,(C)ey„+(C;fc)ey -(C;fc)e 2 ^ (<:/e) d( 

= ^ [ +1 d c (qj(Oev 0+ (C; k)e Vo .(C; k))^-e 2 ^^ d( + V{ 1 

= ^ £ +1 a 2 fe(C)ev + (C; fc)ey _(C; k)) [^j ' e 2 ^^' dQ + b{> 1 + 

with the following boundary terms 

b i = WT] fe(«r + i)eyo+(ai+i;fc)eyo-(«/+i;fc))e 2 ^' (a!+l/e) 
- qj (at)e Vo+ (ai; k)e Vo -{a V , k))e 21 ^'^) , 

- d c (q 3 (C)e Vo+ (C;k)e v ^(C;k))\ c=< e 2 ^'^) . 

Now, one has 

d 2 q 

d l{qj{x)e Vo+ (x; k)e Va -{x; k)) = -j-£(x)e Vo +(x; k)e Vo -(x; k) 

+ 2^-(x)d x (e Va+ (x]k)e Va -{x;k)) + 2q j (x)d x e Vo +{,x\ k)d x e Va -(x; k) 
+ qj(x)((d 2 ev +(x;k))e Vo -(x;k) + e Vo+ (x;k)d 2 e Vo -(x;k)). 

The first three terms are piecewise-C 1 , so that oscillatory integrals predict that 

l {^^ ev ° + ^ fc ) e v -(C; k) + -^(()d c (e Vo+ ((:; fc)ev -(C; k)) 

+2 qj (Od ( e Vo+ (C; k)d ( e Vo -((; fc))e 2 ^' (c/e) d( = 0(e). (5.28) 



SCATTERING FOR OSCILLATORY POTENTIALS WITH STRONG SINGULARITIES 33 



For the fourth term, we use the fact that ey + and ey _ satisfy ( — g^j + Vb — fc 2 )u = 0, 
so that one has, with flj = {xq, • ■ • , xn-i} fl (oj, Oj+i), 

Z € 2 2 ' 

*l [&] = 8fc?r 2 -2 E 2c i*( a:: e Vb+(^; fc ) e Vo-( a; »; fc ) e ^ + b{ 1 + bl' 1 + O (e 3 / f) 

o 



7T 2 7 2 



^ [d x ( qj (x)e Vo+ (x; k)e Vo .(x; k))} + b{ 1 + O (e 3 /j 2 ) . 



Finally, we have ti [q e ] = X^=o J2\j\>i *i* [?e] + C(e 3 ), and one recovers immedi- 
ately terms of the expansion of Theorem 12.11 



M-l 



£ £ ft* = et\ and 

'=0 |j|>l 
M-l . 2 

E E^I-I E [dM^v a + {x-k)e Vo ^x-M] a e 2 ^^ 

2=0 \j\>l XiEilj ° J 

so that ti[g e ] = e*f (fc) + e 2 t|(fc) + C(e 3 ). 

The small e asymptotics of t2[q e ,q e ]: 

Let us assume that £ is fixed outside supp T4j n£ , , and outside the discontinuities 
of qj,d x qj (this particular case arises for a finite number of values of £, and there- 
fore brings no contribution to the transmission coefficient, when integrated). Then 
integrating by part leads to the following expansion for e small: 

7/(0 = -e Vo+ (C,k) t d z {q 3 (z)e Vo+ {z-k)e v ^{z-k)) -±-^'h) dz 
J —oo Hit j 

E \ ey 0+ (C;fc) f - [ C a 2 ( (Z ,(z)ey 0+ (z;fc)ev^(z;fc)) e 2 ^(*A) dz 



4ir 2 j 2 
+ 



d z (q 3 (z)e Vo+ (z; k)e Vo .(z; k)) e ***iU'j\ C 



The first term, treated as previously and using the fact that the functions qj, ey 0+ , 
and ey _ arc piccewise-C 3 , brings a contribution of order 0(e 3 ). 

Now, using the same analysis on and the Wronskian identity (|3.8[) . one 

obtains the following expansion for the — d( integrand of f|5 . 22|) : 



I J i(0+im = (ey 0+ (C;fc)a c fe(C)ey 0+ (C;fc)ey _(C;fc)) 

- ey _(C;fc)a C fe(C)ey 0+ (C;fc)ey -(C;fc))) + 0{e 3 ) 
= -2ikt%° m qj (0 e Vo+ (C;k) e 2 ™^ + 0{e 3 ). 
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Therefore, one has 

1 1 f 2 I + °° 

|j|>1 J 00 |m|>l 

(-2»fctr%(C)ey + (C;fc)e 2 ^ (f/e) ) dC + ^(e 3 ) 
= £2 8^ / + °° £ g - j( ^ (C) ev -(C;fc)e yo+ (C;fc) d( + 0(e 3 ). (5.29) 

One recovers finally: t 2 [&, q e ] = e 2 t>> om + 0(e 3 ). 
Estimate of t e rem : 

Using Proposition I5.3I with Theorem I5.5I yields: 

Proposition 5.6. Let K denote a compact subset of K, satisfying Hypothe- 
sis (G). Introduce for k G K 

tl em (k) = t'(k) - t h Q ° m (k) - et\(k) e 2 ( t h 2 ° m (k) + t%{k) ) . (5.30) 

Then we have 

1. IfV has compact support, thent% em (k) = 0(e 3 ). 

2. IfV is exponentially decreasing, then t^. em (k) = 0(e 3 ~). 

3. If(x) p V G L 1 , p > 9, then there exists 2 < /3 < 3 such that t e rem (k) = £>(e /3 ). 
The proof of Theorem 1 2. II is now complete. 

Appendix A. The numerical computations. In this section we outline the 
numerical method we used to obtain results displayed in figures 12.11 and 12.21 

We approach the computation of t(k), the transmission coefficient associated with 
the potential V(x), by numerical approximation of the function 

u(x;k) = —p-ev-ixik), 
t{k) 

where ey_(i;fc) denotes the distorted plane wave generate by an incoming wave from 
positive infinity; see p. 71) . We rewrite the equation 

d 2 \ 
-— + V(x)-k 2 \u(x;k) = 0, 

equivalently in terms of the variable U(x; k) = (u(x; k), d x u(x; k)) T as the first order 
system 

Note that if V is assumed to have compact support (supp V C [— M, M] with M > 0), 
then 

U(x;k) ee ( e ik l-ik*) iotx<-M, (A.2) 

/ r r (k) n jkx j 1 „-ikx 



U(x; k) ee | ^fg + I for x > M. (A.3) 



t(fe) e t(k) e 
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Starting with the initial data given by (|A.2I) , wc numerically solve the system of first 
order ODEs defined by (|A.1[) up to x > M, and (|A.3[) allows to recover the desired 
value of t(k). At the location of the singularities x = Xj, the jump conditions (13.31) 
allow to obtain U(x+; k) from U(x—\ k) via a transfer matrix. Between the singular- 
ities, one approximatively solves (|A.1[) using for example Rungc-Kutta formulae. We 
used the Matlab solver ode45; see [14] for more information about the Matlab ODE 
Suite. 

We conclude this section by stating the precise functions and parameters used to 
obtain the plots displayed in figures 12.11 and 12.21 

For the case when Vq has singularities, as in the left and center panels of figure |2~T1 
we set 

V = V sing {x) ee 40 (6{x) + 5{x - 0.5) + S(x - 1)) . 

Otherwise, we set 

V = V reg (x) ee 40 (S p (x) + S p (x - 0.5) + S p (x - 1)) , 

with 5 p (x) ee -^=e~ x ' ' l p " the smoothed out approximation. One has p = 0.1 for the 
right panels of figures |2~T1 and l2~2l and respectively p = 0.01 and p = 0.001 for the 
center and left panels of figure l2~2l 

We set q e {x) = f{x)sin{2-Kx/e), with f(x) ee for x e R \ [-2/3; 2/3], and 
elsewhere 

!f{x) = 40 in the discontinuous case (left panel of figure [2~Tj) , or 
f(x) = 40e u-2/3)(x+2/3> m the smooth cases (all other panels). 

Finally, we set k = 5.5, since it corresponds to a case where tQ° m (k) approaches unity 
when Vq = V sing . 

Appendix B. The Jost solutions. In this section, we provide a construction of 
the Jost solutions and a rigorous derivation of their properties, including bounds that 
are used in the proof of Proposition l5.4[ Appendix [Cl We recall that by Definition [321 
the Jost solutions are the unique solutions f±{x; k) of 

(H w - k 2 )u ee (^--^ + W(x) -k 2 ^j u = Q. (B.l) 

such that f±(x:k) = e ±lkx m±(x;k) and 

lira m±(x;k) = 1. 

x— > ±oo 

The existence of Jost solutions for regular potentials W € L 1,3 / 2+ (]R) is established 
in [3]. The generalization to potentials allowing a singular component 

W = Wreg + W s i ng , with 

W reg € L 1 ' 3 / 2 +(M), 

JV-l 

Wsing = Cj S(x — Xj), where Cj,Xj <S K, Xj < Xj+\. 

3=0 
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can be found in [5]. 

As an intermediate step of the proof, one introduces an equivalent definition of 
the Jost solution, as solutions of integral equations. In the case where W is regular, 
one has 



m+(x;k) = l + D k ((-x)W(Om+((;k)d(, (B.2) 

J X 

m_(x;k) = l + f D k (x-OW(Om-((;k)d(, D k (x) = f ' e 2i ^dQ. 

J -oo JO 

If W has regular and singular components, we work with a variant of equations (IB.2|) : 
D k (( - x)W{()m + ((;k)d( + ^ D k (xj - x)cjm + (xj; k), 

Xj>X 

m_(x;fc) = l+ / D k (x - ()W(Qm-(C,k)d( + ^ D k (xj - x)cjm + (xj-,k). 

"'-° Xj<X 

From these integral equations, one deduces 

\m + (x; k)-l\< 1 + 1 "+^ X ' 0) + \»\W(s)\da, 

\m.(x;k)-l\< 1 + 7^p 0) fy + \s\)\W(s)\d S . (B.3) 

Then, since m + satisfies 

d x m+(x;k)= / e Mk{t - x) W{t)m + (t;k), and 

J a; 

d k m+(x;k)= D k {t-x)W{t)d k m+(t-k) + / d k D k (t - x)W(t)m + (t; k), 



one obtains easily the following uniform bounds 

\m + (x;k)\<C{x), \d x m+(x;k)\<C, 

\d k m+(x; k)\ < C(x) 2 , \d x d k m+(x; k)\ < C(x), (B.4) 

where C is independent of k. The same bounds clearly hold for m_(a:; k). 

Appendix C. Proof of Proposition 15.41 This Section is dedicated to the 
proof of Proposition [5]11 namely 

Tr Vq (fc) = (-Do) xRv (k)x (Do) is a bounded operator from L 2 to L 2 . 

This result has been proved by in [3], for Vo = 0, and spatial dimensions n = 1, 2, 3. 
We generalize this result in the one dimensional case for Vo = V reg + V S i ng as in (|2.ip , 
so that singularities in the potential arc allowed. 

Our proof requires the use of the generalized Fourier transform, described in terms 
of the distorted plane waves. We introduce 

mT . n _LJ ev +(x;C) C>0, _ 1 [ t(CW(x;C)e" c C > 0, 
{ '^ } e Vo _(x;-C) C<0, ~ VfcF I t(-C)m_(x;-C)e^ C < 0, 
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with m + (x;C) — > as x —> oo and m_(x;£) — > as x — > — oo. 

Then J 7 and J 7 * the distorted Fourier transform and its adjoint arc defined by 



L? -^L 2 

<p ^7[<t>] (0 = f+™ dx, 

L 2 -> L 2 



7 



One has the following property: 

P c <b = 7*7& 

where P c denotes the spectral projection onto the continuous spectral subspace asso- 
ciated with the operator 

H = -d 2 x + V . (C.l) 

To construct a smoothing operator which commutes with functions of H , it is 
convenient to introduce, using the distorted plane wave spectral representation of H : 

(Dvo) I = (I-A + V ) 1/2 f = [ (r,)7[f](ri) ^(x;rj) d V (C.2) 

JR 

Therefore, one has 

T Rvo = (D ) (Dy,)" 1 (D Vo ) x Rv„ (k) X (D Vo ) (Dy,)" 1 (D ) (C.3) 
= (D ) (D^y 1 o f Rvo o (Dy,)" 1 (Do) . (C.4) 

There are thus three terms to estimate. In order to deal with (Do) (Dy ) 1 and 
(Dy ) -1 (Do), we introduce the classical wave operator, W and its adjoint W* , defined 

by 

W=a - lim e ltH e- ltHo , (C.5) 

t—>oo 

W* = s - lim e UHo e- ttH P c , (C.6) 

t— yoo 

with H = —d 2 + Vq and H = —d 2 . The wave operators have the property to 
intertwine between the continuous part of H and Hq, so that for any Borcl function 
/: 

f(H)P c = Wf(H )W*. 
Especially, one has (Dy ) = W (Do) W* , so that 

(D )(D Vo y 1 = (D )W(D y 1 W*. (C.7) 

Let us state the following result, that has been introduced in [?] and extended 
in [5] to potentials Vo = V reg + V S i ng as in (|2.1[) . thus allowing Dirac delta functions: 

Lemma C.l. W and W* have extensions to bounded operators on H k , for k 
I . i i. 1. Using this last result and the known fact that (Do) 5 is bounded from H k 
to H k ~ s , we obtain directly from (|C.7j) that 

(Dq) (Dy ) 1 is bounded from L 2 to L 2 . 
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Similarly, 

(ZVo) 1 (Dq) is bounded from L 2 to L 2 . 

In order to deal with the last term of (|C.4[) . we decompose Tr v as a sum of four 
operators, commuting (-Do)- 

Tr Vo = (D Vo ) xRv (k)\ (D Vo ) 

= (x (D Vo ) + [(D Vo ) ,X]) RvM ((D Vo ) X + [X, (Dv )\) 
=X (D Vo ) R Vo (k) (D Vo ) X + ([(D Vo ) , X ])(Rv (*) (D Vo ) x) 

+ (x (D Vo ) R Vo (k))([x, (D Vo )]) + (i(Dv ),x])(Rv (kMlx, (D Vo )}) 
=A I + A<$+A?}+A JII . 

Each of these terms is proved to be bounded from L 2 to L 2 . We treat each term 
separately in Propositions IC.21 IC.41 IC.51 and IC.6I 

Proposition C.2. Aj = x (Dv„) Rv„(k) (D Vo ) x ^ bounded L 2 — >• L 2 , i.e. 

\\ X (D Vo )Rv (k)(D Vo )xg\\ L 2 < C\\g\\ L 2, g E L 2 (M) (C.8) 

First we commute (Dy ) and Rv a - It is obvious that Rq and (-Do) commute, so 
that using the wave operators introduced above (so that (-Dy ) = W (-Do) W* and 
R Va {k) = WR Q (k)W*, with W unitary). 

Ai = x(Dv )Rv (k)(D Vo )x 

= X W (D )W*WR (k)W*W (D )W*x 

= xWRo(k)(D ) 2 W* X 

= xRv„{k) (D Vo ) 2 x- 

Then, applying the identity (Z?y ) 2 = I — A + Vo, one obtains 

Aj = (1 + k 2 )xR Vo (k)x + X 2 - 
Finally, using p. lip together with (|B.4|) . one has the pointwise bound 

\R Vo (x,y;k)\ < C(x)(y) 
with C uniform in k. It follows that for f E L 2 , 



\xRvo(k)xf\ L 2 = 



X {x) / i?y (x,C;fc) X (C)/(0 dC 



c 



< C\ X (x)(x)\ 2 Li If 



L- 



L- ' 



so that Ai is bounded from L 2 to L 2 , with 



\M\l^l> ^ C (\x(x)(x)\lz + kLoo) ■ (C9) 



Before carrying on with estimating the term A^j , let us state the following 



Lemma. 
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Lemma C.3. Let K be defined for (£,7?) G K x K fey 



= ((f) - ( ?? )) y *(C;0 *(C;»7) x(C) 



TTiera K(£,r]) satisfies the following upper bounds: 



c' 



| d v K^,ri) | + | d f Ktt, v ) | < i - *_ - , 



iwi/i i/ie C x and C' constants depending on the function x with 



c x = c ]T 11(0' dlx + ||(C> 2 x|L + ||(C> 2 x|| L . 

V=o " L < c 

= c IE |(C) J ' +1 ^x rl + ||(C) 3 x\\ L , + ||(C) 3 x| 



(CIO) 

(C.ll) 
(C.12) 



Proof. We consider the case where £ > and 77 > 0. The other cases follow 
similarly. Therefore, one has 

K(tv) =((0-(r?» !«,»?), with 



l(^v) = J *(C;0*(C;»?)x(0 ^C = y e^-^(0m + (C;0^)™+(C;7?)x(C) C 

(C.13) 

Throughout the proof, we will use extensively the uniform bounds on m + given 
in (|R4]t . 

First, by the uniform boundedncss of and (C) _1 £) m C ancl £j one nas 

|^,r?)|<|t(OII^)l / KcrWC;??) (crWOOKO'MOK <c|I(C) 2 xIIl- 

(C.14) 

For |r? — £| > 1 we write 

i&n) = ( . fa * e))2 *W?) / ( ^ e iC( "-° ) ^7(C^) m+(C;C) x(C) dC 

(C.15) 

= (i(r,-i;))* *®t(v)Je«w ^ (m+iCv) m+(C;0 x(C)) dC. (C16) 



The most singular terms in the integrand of (|C. 16|) are those containing 3j?m + . In 
particular, recall the relation 9^m+ = — 2ikd x m + + Vom + , where Vq contains Dirac 
mass singularities. Thus, for |£ — 77 1 > 1 we have 



m,v)\ < c (E ||(o j ^x 



1(0 x 



1^. + ll(0 2 x|| io 



(C.17) 
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Applying (jCll) for \r) - f | < 1 and (fCTF)) for \r] - f | > 1 yields 

|/«,»7)l < C x 



1 + 

Finally, since 77) | = 77) | - (rj)\ < C |J(£, 7?)| |£ - r)\, multiplication by 

\C-T}\ implies (jCTTT]) . 

Using the same method as previously, one obtains similarly 

1 



\d v m,v)\ < C' 



1 + \t ~ V\ 



Finally, one has \d v K(^,r])\ < \d v I(^,r])\ — (ri)\ + J7)|, so that we deduce 
the first part of (|C12|) . By symmetry, one obtains the same estimate for d^K(^,rj), 
which concludes the proof of Lemma IC.3I □ 

Proposition C.4. Ay/ = [(D Vo ) ,x\Ry ( D v ) X is bounded L 2 -> L 2 , i.e. 

\\{(D Va ),x}Rv (Dv )x9\\ L z < C\\g\\ L 2, g g L 2 (M). (C.18) 



Proof. Our strategy is as follows. We view the operator A// as a composition of 
two operators 

4/ = l(Dv ),x) ° Rv (Dv )x 

and first find a representation of each operator with respect to the distorted Fourier 
basis. We then directly prove the boundedness of A// : L? H> L 2 using this spectral 
representation and an appropriate frequency localization argument. 

In terms of the distorted Fourier transform, one has 
l(D Vo ) ,x]f(x) = [{D Vo ), X ] (J ^(x;vmf](v) dqj 

= J nfKv)((Dvo)(x(xMx;v))-x(Dv a )*(x;vi) di], (C.19) 
Now, since (Dy ) ^(x;rj) = (r/)4 , (x; rf), one has 

(D Vo ) (x(x)*(x; V )) = J *(x; j (D Vo ) (x*(s v))(C)WiO d( <% 
= J J x(0* (C; rj) (Dv )W®(0 d( 
= J J X (C)*(C; V)(0W^) d( «■ 

Therefore, we finally deduce 
[(AO , x]/(z) = I .F[/]fa) ^ *(Ci ^7)^(C0x(C) «0 - fo» rfC ) 

= / f ({0-{v)) f*M*((;i)x(CW/]W«. (C20) 

J? J r) JC 
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To represent the operator Ry a (Dy ) x m terms of the distorted Fourier basis we 
note: 

F[R Vo (Dvo) xg}(v) = / *(*;»?) (Rvo (D Va ) xg){z) dz 

Jr] 

= / (R Va (D Vo )^{z;r ] ))x(z)g(z) dz 

J Z 

(v) 



i] 2 — k 2 
(V) 



V(z;r))x{z)g(z) dz 



F[xg](v)- (C21) 



rj 2 — k 2 

Combining (|C.20|) and (|C.21|) . one has 

l(D Vo ), x ]Rv (D Vo )xg(x) = f f «0 - (v)) [ WiO^(C,v)x(0 d( 

J£ Jr) JQ 

-fiL^F{ X g}(v) d v dt, 

*(x;0T u [g}(Z)dZ. (C.22) 



By the Plancherel Theorem, the L 2 estimate of A^j is equivalent to the bound 



T n [g] 



\L' 1 



r]J( 



*(C;0 x(0 dC (<€) 2 ^ <7?) rbsKv) dv 



<C\\g\\v 



(C.23) 

We now proceed with a proof of (|C.23|) . First we define tp\ K \<$ to be the positive 
smooth function satisfying 

l P\k\<8 equal to one for |k| < <5o/2, zero for |k| > Sq and symmetric about k = 0. 

(C.24) 

We use tp to localize at frequencies near 77 = ±fc and frequencies away from rj = ±fc. 

T u [g] = T^M + Tj-lAg] (c.25) 

where 

T}L[g}(0 = [K&V)-rKs [^-nr,\-\k\\<s (v)}nxg}(v)dv, (C.26) 

Jr] T) — K 

TneaMiO = I Kfov) v M -W<*M ^MM *n, (C.27) 

Jr] V K 

with K defined as in (fCTO) by 



I 



K&r)) = «£>-<??» y *(C;0*(C;»7)x(C)dC- 

Bound on Tf* r \g](£): We bound the expression 

TjlAgm = Jk^t,)-^^ [l-^ lvHkll<So ( V )]nxg}(v)dv. (C.28) 
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By Lemma [C. 31 K(£, rf) satisfies the following pointwise bound, which is valid for 
all £, r) 6 R: 

\mrj)\< C x ^— -1— j.. 
Recall now the special case of Young's inequality: 

\\h*gh < \\ h h hh ■ 

This, together with the pointwise bound of K(^,r]), yields: 



T "«M II 2 



<■'/> 



I] 2 — k 2 



[1 - f\\ v \-\k\\<s (v)] I F[X9](v) I dr l 



l 

i 



n 

2 



£?. 



{''/) 



£1 



< c x \\x\\l°° UWl 2 - 



Bound onT4 7 ar M(0 



rpII,E 

near 



^'^) ^2 _ fc 2 ^e<lk|-|fc||<^o(^) -^[XSK 7 ?) 



1 [ jrit m \ < \ f ^lxgKv) nxg](v) \ , 

2fc y/^^^^iN-WK*,^) (^-^ ^^J^ 



where 



A £ (^) = ^ fa) Ve< M -\k\\<s (v) K(t,ri). (C.29) 
and rj) is displayed in (|C.10[) . Note that by Lemma [C. 31 

|A e (£,77)| < C X l + ^_ Ve < M . M<s M- ( C - 3 °) 

We bound the first term in the above expansion of T^^ r . The second term is 
treated similarly. We have 

J A £ (£, ?7 ) = ^(fl + + where 



\z/t i,w ^IxgKv) * , 

„_fc 1 l'7-fc|>«o/4 









/ (A e &7?) 




/ A e (^,r,) 
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One bounds Tl e using (fC30)) by 



TV 



< 



4C, 
4C* V 



We<\\ v \-\k\\<s a (r])T[xg]{r])\\L^ 



— llv £ <n»7i-ifcii«solU=llxlU°° HsIU 2 - 



(C.31) 



Moreover, we have 
A^r^-A^fc) 



77 — fc 



< 



9»j A e (^,7?) |^=^ e { e <|^_jfe|<5 /4} 

< l|,_ fc |<i ( I 9, ^,77) I + (C.32) 
From the estimates of Lemma |C.3[ and using Young's inequality, one deduces 

II £ e IU« < c£ ||xlU-NU». (C33) 

We treat the singular integral S £ as follows. By antisymmetry of the function 
(r) - k)- 1 l £ <\ n -k\<s /i(v) w e have 



/ 



1 



L e<l'?-fc|<<5o/4 



T) — k 

Moreover, we have 
77 — k 



F\X9](v) dri 



L e<|r;-fc|<<5o/4 



F[ X9]{ri) ~ F[X9\{k) 
77 — k 



dr\. 



< 



dri ^IXdliV) \T l =f)€{e<\i 1 -k\<Sa/'l} 



By the uniform boundedness of (() 2 d v m + ((,ri) and (C) 2 d v m-(( 7 r]) in 
have that 



(C.34) 



x K, we 



\d v (F\x9Kv))\ 



9 v nC;v)x(0g(0 d( 



Therefore, 



< sup |(C)- 2 ^*(C;»7)IIKC> 2 xIIl?NU» < c|Kc) a xlU?NI 



< c \\(Q 2 x \\ Ll \\g\\ L 2, 



L- ■ 



1 e<\ v -k\<8 /4(v) 7 FixgKv) dr) 



from which it follows that 

1^(0 I < C \\(Q* X \\ L ,\\g\\ L , |A £ (£,fc) I < 



\\gh». (C35) 



x l + l£-fcl 

Thus we have from (ICT3lT) . (I(l33l) and (jH^Hjl . 

|| Tii ar [g] || 2 < C' x \\g\\ 2 . 

Using the estimates of Tj^ r [g] and T™ ar [g] yields (^23]) . Therefore, is 
bounded from L 2 to L 2 . This completes the proof of Proposition IC .51 □ 
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Proposition C.5. A}} = \ (Dv ) R Vo (k))([x, (D Vo )}) is bounded from L 2 to L 2 . 
This follows from Proposition IC.4I and duality. 

Finally, we consider the operator Am = [(Dy ) ,x] ° Rv (k) ° [x> (Dv )]- 
PROPOSITION C.6. The operator Am is bounded from L 2 to L 2 . 
Proof. By (|C.19|) . one has 

A ln [g](x) = f *(x;0 / ((0 - (r,)) f W^)*(C;vM0 d( 

--L^F[[{D Vo ),x\g](v) dvd( 
= f *(a:;0 / K^,r,)^—T[[(D Vo ), X }g](v)dr,d^ 

J£ Jr\ J] — K 

= f *(x;0 f Kforj) ^-rr 2 f K( V ,0)F[g}(0) d9 dr, d£ (C.36) 

J£ J 7) 1] — K Jg 



By the Plancherel Theorem, the L 2 estimate of Am is equivalent to the bound 

1 



T ul [g] 



2 -k 2 



K(r),6)T[g}(0) d6 dr] 



< C \\g\\ L 2 



(C.37) 

We now proceed with a proof of (|C.37|) . We use (fi\ K \<s , defined as in (|C.24|) . to 
localize at frequencies near rj = ±fc and frequencies away from r\ = ±fc. 



(C.38) 



where 



TjtMit) = m,ri) (l-V|N-|fc||<*o(»/)) 



rf- — k 2 



T near[g}(0 = / K(£, Tj) lf\ |,|_ | fc | | <{ „ (r?) 



1 



r) 2 — k 2 



K(r],0)T[g}{0) d9 drj 

(C.39) 

K^q.e)^^) d9 dr) (C.40) 



Bound on Tj^[g]{^): We recall Lemma [C.3[ stating that K(^,rj) satisfies the fol- 
lowing upper bound: 

I I < c XY ^— y 

with C X = C fo =0 \\(CV 8lx\\ Ll + || (O 2 x|| L i + ||(C> 2 x\\ L Sj- Therefore, one 
has the pointwise estimate 



K(r,,6)F{g}(0) d9 





1 






1 


< 




T 2 


\\\v< 




1 + \r) - 0\ 


l + l-l 



Iff Ilia 
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Moreover, for ||r?| - |fc|| > S , one has (l - <P\\ v \-\k\\<8 (v)) \v 2 ~ k 2 \ 1 G L 1 . 
Therefore, by Young's inequality, 



ra<?]|| L2 <c 



1 


2 


l+l-l 


L 2 



I 1 ~ fWvl-MKSoiv)) 



1 



|r/ 2 - fc 2 | 



IIsIIl* <CxlMI £ « 



Li 



Bound on [«?](£) 

r^; 8 [s](0 = 



i 



K(£,r)) <P\\ v \-\k\\<6o(v) _ 2 p 
A e (^,r/) — ^- / tf(»7,0).F[ff](0) dfofej 



/<(?7, 6')J r [ ff ](6») d<9 dry 



r\ — k 



A £ (£,r/) 



1 



rj + k 



K(r,,e)T[g}(0) d8d v , 



withA e (^,?7) = ^ A'(C,r?) ^ E <|M-|fc||«j (T/)- 
Note that by Lemma [C .31 

I A e (£,r,) | < C x 



1 + 



t Ps<\\ri\-\k\\<S (v)- 



(C.41) 



We bound the first term in the above expansion of T^^f. The second term is 
treated similarly. We have 

A e (^) = + + ^(0, where 



S £ (0 =A e (£,fc) 



1 



V - k J 



n%0 = / A e (t,v) 



TJ — k 
1 



K{il,0)T[g\{e) d.6 l e <| TJ _ fe |< 5o/4 d-q, 
K (v,9)^[9}(G) d9 l e <|, ; - fe |<5 /4 dr], 



7} — k 



K( Vl 6)T[g](d) d9 l W _ fe |> 5o/4 dr,. 



As in the proof of Proposition IC.41 the kernel of the integral operators defining 
£ e and 1Z 6 are non-singular, and we have uniformly in e: 



\\£ £ \\l> + \\K s \\ l , < C'\\g\\ 



L- ■ 



We treat the singular integral S £ as follows. By antisymmetry of the function 
(rj - A:) -1 l e < | »7— fc| <<S" /4 (»7) w c have 

s(0 = A(e,fc) f i e <\r,-k\<6 0/ M I HgW) de dr, 



T[g](6) d9 drj. 



1 £<|7 ? -fc|<,So/4( 7 ?) 

Moreover, Lemma [C . 3 1 leads to 

K( v ,e)-K(k,t 



r] — k 
K(r,,e)-K(k, 
r\ — k 



tj — k 



< 



d v K(t],6) l 7)= jjg{e<| 7; -fe|< < 5 /4} 



: CL 



1 



x i+i?j-^r 

(C.42) 
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Therefore, by Cauchy-Schwarz inequality, 
K(r], 9) - K(k 



i] — k 



F[g]{9) dO 



< C'\\g\\ 



L 2 ! 



from which it follows that 



\s%0\ < |A e (€,fc)| ||i e <|,_ fc |< io/4 ||Li \\g\\ L * < c x 1 + | e _ fc i WaWv- 



(C.43) 



Thus we have || T^ r [g] \\ 2 < C' x ||g|| 2 . Using the estimates of T}£\g] and 
^nearidl yields (|C.37|) . Therefore, Am is bounded from L 2 to L 2 . This completes the 
proof of Proposition [CT6l and hence the proof of Proposition [5T4] □ 
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